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ABSTRACT 

We present new measurements of stellar velocity dispersions, using spectra ob- 
tained with the Keck Echellette Spectrograph and Imager (ESI) and the Magellan 
Echellette (MagE), for 76 Seyfert 1 galaxies from the recent catalogue of Greene 
& Ho. These objects were selected from the Sloan Digital Sky Survey (SDSS) 
to have estimated black hole (BH) masses below 2 x 10^ Mq. Combining our 
results with previous ESI observations of similar objects, we obtain an expanded 
sample of 93 galaxies and examine the relation between BH mass and velocity 
dispersion (the Mbh — c"* relation) for active galaxies with low BH masses. The 
low-mass active galaxies tend to follow the extrapolation of the Mbh — o"* relation 
of inactive galaxies. Including results for active galaxies of higher BH mass from 
the literature, we find a zero point a = 7.68 ± 0.08 and slope of /3 = 3.32 ± 0.22 
for the Mbh — cr* relation [in the form log Mbh = a + /31og(cr*/200 km s~^)], 
with intrinsic scatter of 0.46 ± 0.03 dex. This result is consistent, within the 
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uncertainties, with the slope of the Mbh — cr* relation for reverberation-mapped 
active galaxies with BH masses from 10^ to 10^ Mq. For the subset of our sam- 
ple having morphological information from Hubble Space Telescope images, we 
examine the slope of the Mbh "~ cr* relation separately for subsamples of barred 
and unbarred host galaxies, and find no significant evidence for a difference in 
slope. We do find a mild offset between low-inclination and high-inclination disk 
galaxies, such that more highly inclined galaxies tend to have larger a^, at a given 
value of BH mass, presumably due to the contribution of disk rotation within 
the spectroscopic aperture. We also find that the velocity dispersion of the ion- 
ized gas, measured from narrow emission lines including [N II] A6583 A, [S II] 
AA6716, 6731 AA, and the core of [O III] A5007 A (with the blue-shifted wing re- 
moved), trace the stellar velocity dispersion well for this large sample of low- mass 
Seyfert 1 galaxies. 

Subject headings: galaxies: active — galaxies: dwarf — galaxies: nuclei — galax- 
ies: Seyfert 



Introduction 



Dynamical studies of local galaxies over the past decade have firmly established that 
supermassive black holes (BHs) are present in most (possibly all) galaxies with massive 



bulges, and that the black hole mass tightly c orrelates with 



locity dispersion (jFerrarese fc Merritt 



2000 



j ulge mass and stellar ve- 



Gebhardt et al. 



2000 



Merritt fc Ferrarese 



2001; Tremaine et al. 2002 



Marconi fc Hunt 



2003 



Haring fc Rix 



2004 



Giiltekin et al. 



20091 



These correlations suggest coeval growth of the galaxy bulge and the central BH. The 
-Mbh — 0"* relation is import ant both as a fundamental b enchmark against which galaxy 
evolution models are tested ( iHaehnelt fc KauffmannI 120001). and as a key input to calcu- 



lation s of the density of black holes in the Universe ( iMarconi et al.l 120041 : IVolonteri et al. 
2008h . 



The low-mass end of the mass function locally is where models of primordial BH seed 
formation show pronounced differences, thus improving the observational constraints on the 
low-mass end of the Mbh — o"* relation can provide important constraints on the models. 
Theoretical work has prop osed two types of seed models: "light seeds" (IMadau fc Rees 



2001 



Volonteri et al. 



20031) formed as remnants of Po pulation HI stars, and "heavy seeds" 



( iKoushiappas et al.l 120041 : iLodato &: NatarajanI 120061 ) from direct collapse of massive gas 
clouds in primordial halos. These different classes of models should result in different de- 
mographics for BHs in low-mass galaxies at the present epoch. Low-mass seeds lead to a 
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wide scatter in BH masses in low-dispersion galaxies (with some galaxies hosting very low- 
mass BHs) but a relatively high occupation fraction of BHs in low-mass galaxies, while in a 
heavy seed scenario, the BH occupation fraction is low but the minimum BH mass is larger 
( iVolonteri et al.ll2008t IVolonteri fc Natarajanll2009l ). Unfortunately, observations have not 
yet been sufficiently sensitive to distinguish between these scenarios, and additional obser- 
vations to constrain the Mbh — c* relation at low masses are needed. Empirically, the slope 
and scatter of th e Mrh — cr* re lation are still subject to debate, particularly at the low and 



high mass ends (jWyithd 120061 ) . There is some recent work that suggests that the relation 
may be different for different types of h ost galaxies, i.e., for barred ys. unbarred galaxies 



and for classical bu lges vs. pseudo-bulges ( lHij|2008l : iGraham &: Lill2009l : iGiiltekin et al.ll2009 



Greene et al.ll2010[ ). These differences may be more pronounced at lower masses and cause 
scatter in the overall Mbh ~ cr* relation. However, due to the difficulty of obtaining direct 
stellar-dynamical measurements of BH masses in low-mass galaxies, much of the information 
on BH demographics at low mass comes from active galactic nuclei (AGNs). Recent Hubble 
Space Telescope (HS T) imaging of low-mass Seyfert 1 galaxies selected from the SDSS by 
Greene fc Hoi (120041 ) also shows some evidence for a change in the Mbh — -^^buige relation at 
low mass, probably indicative of a different mode of BH growth in th ese objects compared 
with higher- mass galaxies having classical bulges ( iGreene et al.l 120081 ). In this paper, our 
goal is to expand the sample of low-mass AGNs having both black hole mass estimates from 
single-epoch spectroscopy and direct measurements of stellar velocity dispersion, in order 
to investigate the low-mass end of the Mbh — cr* relation in more detail than was possible 
previously. 



Greene fc Hoi ( 12004 . l2007d . hereafter the GH07 sample) searched the SDSS database 
and presented a large sample of AGNs with low-mass BHs (Mbh < 2 x lO^M©; see also 



Dong et al.ll2007l . for related work). The for the GH07 sample could not be determined 
from the SDSS spectra, because the instrumental resolution of ~ 70 km s~^ sets a practical 
lower limit to the that could be measured, and also because the S/N of the SDSS spectra 
was insufficient or the continuum was dominated by AGN emission rather than by starlight. 
Here we present new measurements of cr* for objects in the GH07 sample. These new 
measurements provide a useful way to examine BH demographics, even if the individual 
BH masses are not highly accurate. In addition, our work provides a very large sample 
to examine relationships between gas and stellar kinematics for nearby AGNs. We present 
the sample properties, observations, and data reduction in Section |2l and measurements of 
(T* and line-width in Section |3l We discuss the Mbh — o"* relation in Section H] and the 
narrow-line properties Section El with conclusions in Section El 
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Observations and Data Reduction 



The objects o bserved at Ke c k and Magellan were selected from the sample of low-mass 
BHs presented by iGreene fc Hoi (l2007d ). They selected 229 broad-line active galaxies with 
Mbh < 2 X IO^Mq from SDSS DR4. They estimated the BH masses with the single-epoch 
virial method, using the Full Width Half Maximum (FW HM) and lumi i iosity o f the broad 
Ha emission line, following the methods first described by lGreene fc Hoi (j2005bl ). The virial 
calculation det ermines the r adius of the broad-line region (BLR) using the radius-luminosity 
relationship of iBentz et al.l ( l2006l ). 



The o bjects were observe d during 2008-2009 using the Echellette Spectrograph and Im- 



ager (ESI; 



Sheinis et al.ll2002l ) at the Keck-II telescope and the Magellan Echellette (MagE; 



Marshall et al.ll2008l ) Spectrograph at the Magellan II Clay telescope at Las Campanas Ob- 
servatory. For each observation, the spectrograph slit was oriented at the parallactic angle. 
Flux standards and late-type giant stars (F4-M0) for use as velocity templates were observed 
during each night. Details of the Keck and Magellan observations and reductions are given 
below. We obtained useful data for 65 of 66 objects observed with MagE and for 13 objects 
observed with ESI; two objects were observed with both ESI and MagE. We also co mbine 
our new sample with the 17 similar objects previously presented by lBarth et al.l (120051 . Here- 
after BGH05). Including the BGH05 sample, our total sample consists of 93 low-mass SDSS 
AGNs. 

ESI - Observations with ESI at Keck were made during the nights of 2008 March 1-2 
UT. We used a 0'.'75 slit width, resulting in an instrumental dispersion of o"j ~ 22 km s~^. 
The spectra cover the wavelength range 3800-10900 A across 10 echelle orders, and the 
dispersion is a constant 11.5 km s~^ pixel"^ in velocity. The exposure times for individual 
objects ranged from 900s to 3000s. One-dimensional spectra were extracted within a 1" 
extraction width, and wavelength- and flux-calibrated, with correction for telluric absorption 
bands, using standard techniques following the same methods we have previously used for 
ESI data (BGH05). 

MagE - Observations with MagE at the Magellan II Clay telescope were carried out 
on the nights of 2008 April 9-13, 2008 August 28-30, and 2009 January 25-27 UT, using 
a 1" slit width, giving an instrumental dispersion of ~ 26 km s~^, as measured from the 
arc lamp spectra. The spectral coverage is approximately 3200-10000 A across 15 echelle 
orders, with a nearly constant dispersion of 23 km s~^ pixel"^. Exposure times ranged from 
1800-7200 s and were typically 5400 s. One-dimensional spectra were extracted from the 
CCD images and wavelength-calibrated using the MAGE_REDUCE package kindly provided by 
George Becker. The extraction width is typically 1'.'5 to 3". T he MAGE _ R,EDUC E package 
uses techniques for rectification and sky subtraction developed by iKelsonI (120031 ) . Optimal 
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extractions (lHorndll986l ) were used for exposures of the bright stars, while for the galaxies 
a simple boxcar extraction was used in order to avoid spurious clipping of emission lines 
that can sometimes occur with optimal extractions. Flux calibration and telluric absorption 
correction were applied using the same methods used for the ESI data. Finally, cosmic rays 
were removed when combining multiple exposures. 

Two objects, SDSS J093147.25+063503.2 and SDSS J131310. 12+051942.1, were ob- 
served with both ESI and MagE, allowing a direct comparison of the spectra; see Figure [U 
To facilitate comparison, the ESI data have been re-binned to the spectral resolution of the 
MagE data, and all the spectra have been re-scaled to have flux density of unity at 5100 
A. The Keck and Magellan spectra appear consistent across the entire available wavelength 
range, which illustrates the consistency of the flux calibration, except for differences in the 
emission lines, which are expected due to the difference in instrumental dispersion of the two 
instruments. For SDSS J093147.25+063503.2, the difference spectrum reveals significant 
residuals only in the peaks of the narrow emission lines, most likely resulting from the dif- 
ference in instrumental dispersion. For SDSS J131310. 12+051942.1, there is also evidence of 
some flux variation between the 2 observations. This could result from the different aperture 
size and different seeing between the two exposures, but could also be in part due to real 
variability of the continuum and broad emission lines. However, the time interval between 
these two observations is 329 days, and on such a short timescale the narrow lines probably 
would not be variable. Since there are residuals in the narrow lines in Figure [T] as well as 
the broad lines, the main cause is probably aperture size and seeing. 



Measurements 



3.1. Stellar Velocity Dispersions 



The stellar velocity d ispersions were measured by a direct fitting method (IBurbidge et al. 



196ll : iRix fc White! Il992l ). in which the spectra of velocity template stars are broadened 
and fitted to the galaxy spectra locally in a specific spectral regi on. A Gaussian p rofile 
was assumed as the line-of-sight velocity distribution. We follow iBarth et al.l (|2002| ) and 
Greene fc Hoi (l2006al ) and express the fitted model spectrum, M{x), as 



M{x) = {[T{x) ® G{x)] + C{x)}P{x), (1) 

where T(x) is the stellar template spectrum, G{x) is the Gaussian broadening function, 
C{x) represents a featureless continuum, and P{x) is a poljTiomial factor. Since our fits are 
performed across a small wavelength region, we adopt a quadratic polynomial for C{x), which 
could also account for some other additive components, such as the "pseudo-continuum" due 
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to Fe II emission originating from BLR of AGN. The low-order multiplicative polynomial P 
allows for the differences between the template and the galaxy, including continuum shape, 
reddening in the galaxy spectrum, and wavelength-dependent flux calibration errors. We use 
a quadratic polynomial for P, since higher-order polynomia ls will tend to fit t he absorption 
features and adversely affect the dispersion measurements (IBarth et al.ll2002l ). Besides the 
velocity dispersion and six parameters for these two polynomials, the redshift of the galaxy 
is also a free parameter. We determined the best-fit parameters by minimizing x^, using the 
Levenberg-Marqu ardt least-squares fitting routine provided by the mpfit package in IDL 
dMarkwardtlboOoh . 



Our collection of stellar templates includes 19 G and K giant stars observed with ESI 
and 22 F, G, K and M giant stars observed with MagE. For each galaxy, we inspect the fitting 
and discard the templates that fail to fit. Then we list the of the remaining templates 
in ascending order, and select the first two-thirds to be well-fit templates. The measured 
velocity dispersion derived with the best-fitting template (minimum x^) is adopted to be 
the best estimate of cr*. Generally, the best fits were obtained with K-giant templates. The 
uncertainty in measurement is calculated as the quadrature sum of the fitting uncertainty 
of the best-fit template and the standard deviation of the measurements of all the selected 
templates. A minimum of six templates are used in each calculation, except for several 
objects, marked by ":" behind the value of a* in Table [H that were not well fit by six or 
more templates. In these cases, all the results of well-fit templates will be accounted in 
estimating the uncertainty. 

The spectral region around the Ca II AA8498, 8542, 8662 triplet (~ 8470 - 8700 A, 
hereafter the CaT region) is ideal for measuring velocity dispersions because the CaT lines 
are strong, not blended with other strong lines, and relatively insensitive to stellar population 
variations. They are also less strongly diluted by AGN continuum contamination than 
stellar features at blue wavelengths. However, for objects with redshift higher than 0.05, the 
CaT absorption features can be affected by both night sky emission residuals and telluric 
absorption bands, making it difficult to obtain useful measurements. Stellar features at bluer 
wavelengths are unaffected by telluric absorption, and we also carried out measurements in 
the Mg 16 region (~ 5050 - 5 430 A) and the "Fe regio n" ( ^ 5250 - 5820 A) r ed ward of Mg Ih 
for me asurements, following iBarth et al.l (120021 ) and iGreene fc Hoi (l2006al ). I Greene fc Ho 
(j2006al ) ran a series of simulations to evaluate the contamination of narrow emission lines, 
including [Fe VII] A5158 A, [Fe VI] A5176 A, and [N I] AA5197,5200 AA, around the Mg Ih 
features, as well as the pseudo-continuum of broad Fe II extending in this region. They found 
that the Fe region containing strong Fe I absorption features resulted in better recovery of 
cr*, since it includes less contamination by coronal emission lines. 
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We tested these two regions by fitting stellar templates to a set of simulated spectra, 
composed as a linear combination of a K2 star, an AO star, and a featureless linear continuum. 
The combined spectra were then broadened by a Gaussian velocity-broadening kernel with 
width ranging from a* = 30 to 100 km s~^ in increments of 10 km s~^. Slight variations in 
the shape of the spectra, representing calibration errors and random errors, were imposed 
on the broadened model spectra, and the models were degraded to S/N = 10, 30, 50, and 
100 per pixel. We also redshifted the model spectra by an arbitrary value comparable to 
the redshifts our observed sample. Then the modeled spectra were fitted using the methods 
described above, and the measured cr* was compared with the input value. We find, for 
6(7^ = [(T(output) — o"(input)]/o"(input), that (5(T*) = 0.02 ± 0.11 from the Mg lb region and 
(5o"*) = 0.03 ± 0.02 from the Fe region for S/N = 30, the median S/N for our observed 
sample. Reassuringly, the results from these test measurements are consistent with input 
values. We find that the scatter in results from the Mg 16 region is larger than that from 
Fe region. The scatter for both fitting regions decreases with increasing S/N, to 0.04 and 
0.01 for S/N = 100, respectively. The larger scatter of the Mg 16 region might be attributed 
to template mismatch. We checked the template spectra and found that for different types 
of stars, the variations in the width of the Mg 16 lines are slightly larger than the width 
variations for the Fe I features. For different type of stars, for example, in the sequence of 
G8 — 7- K3 — 7- K5, the Mg 16 line widths increase, while Fe absorption widths are consistent; 
this may explain at least in part why the dispersion measurements obtained from the Mg 16 
region have larger scatter. 

In practice, when fitting to actual galaxy spectra, we excluded wavelength regions cov- 
ering emission lines including [N I] AA5197, 5200 AA and the high-ionization Fe lines. The 
majority of the objects have several high-ionization Fe lines, including [Fe VII] A5158 A, 
[Fe VI] A5176 A, [Fe VII] A5278 A, [Fe XIV] A5303 A and [Fe VI] A5335 A. We tested model 
fits that excluded and included the Mg 16 absorption features themselves, and found that 
the best-fit template could generally fit the Mg 16 features well for galaxies not dominated 
by AGN Fe emission. This is because most objects in this sample have relatively low cr*, 
so the mismatch of [Mg/ Fe] abundance ratio t hat sometimes affects template fits to ellipti- 



cal galaxies with high a^: (IWorthey et al.lll992l ) is apparently not a significant issue for this 
sample. 

In order to avoid issues of mismatch in flux calibration or spectral resolution across 
echelle orders, we prefer to measure velocity dispersions from a single echelle order (rather 
than from multiple orders that have been "stitched" together), and the specific fitting region 
for each individual galaxy was adjusted to remain within one echelle order. For MagE data, 
we first measured a* from both the Mg 16 region and the Fe region individually. We found 
that if these two regions are in the same order, we obtained highly consistent results from 
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both regions. This is partly due to the hmited wavelength range of each order for MagE data 
and the fact that there is substantial overlap between the Mg 16 and Fe regions. Therefore, 
if these two are on the same order for a given galaxy, we carry out a single fit extending 
across both of these regions, and list the result as o"*(Mg 16) in Table [T] When a separate 
value is listed in the a^(Fe) column of the table, this denotes that the Mg 16 and Fe regions 
fell in adjacent echelle orders and were fitted separately. Figure |2] illustrates some examples 
of both ESI and MagE spectra in the Fe and CaT regions. 

We were able to obtain useful measurements of cr^, for 56 of the 76 newly observed 
galaxies. The two galaxies with both ESI and MagE observations have consistent velocity 
dispersions within the uncertainty measured by two instruments. The remaining 20 galaxies 
are either dominated by AGN emission or have S/N too low to permit a successful fit. We 
also re-measured a* for 15 of the 17 objects from BGH05 (excluding the two objects from 
that sample that were highly AGN-dominated) . The results are consistent with the BGH05 
measurements, within the uncertainties. If we define the deviation between our new mea- 
surements and BGH05 both from the Mg 16 region, as Act* = log cr* (new) — log (T*(BGH05), 
the mean value of Acr^. is less than 0.001 dex and the rms difference is only 0.02 dex. All 
of the (j^, measurements are summarized in Table [U For the galaxies with more than one 
measurement of a* from different fitting regions, the corresponding results are mostly con- 
sistent within the uncertainties. We take the average of the available measurements as the 
best estimate of the stellar velocity dispersion in each object, and list as cr* in Table [H 



3.2. Emission-line Properties 

To decompose the Ha-|-[N II] lines, we first subtracted the underlying continuum. The 
continuum model is essentially the same as that in Equation [1] used to measure a*, but 
fitted over the spectral region surrounding Ha, covering the rest-frame region 6100-7100 
A. For this model spectrum {M{x) in Equation [1]) , the velocity dispersion was constrained 
to lie within three times the la uncertainties around the measured velocity dispersion as 
described in the preceding section. Emission lines were masked out from the calculation of 
for the fits. The continuum fits were carried out on spectra in which the echelle orders 
were "stitched" together, since in some cases the Ha line falls near the end of an echelle 
order. The best-fit continuum model was then subtracted from the spectra to yield a pure 
emission-line spectrum. 



We used multiple-Gaussia n models to fi t the H a and [N II] lines as well as [S II] 
AA6717,6731 AA. We followed iGreene fc Hoi fl2007bf ). using a multi-component Gaussian 
fit to the [S II] doublet to model other blended fines. Up to 4 Gaussians were used to model 
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the [S II] doublet. The velocities and widths of the two [S II] lines were constrained to be 
the same, while the intensity ratio was allowed to vary. For objects with very weak or absent 
[S II] emission, we used the core of the [O III] A5007 A line as the narrow-line model instead. 

The [N II] doublet lines were constrained to have their relative rest-frame wavelengths 
and intensity ratio fixed to their laboratory values of 35.42 A and 2.96, respectively. We 
then fit as many Gaussian components to the broad component of Ha as needed to achieve 
an acceptable fit: starting with a single Gaussian model, new components were added one 
at a time if they resulted in 20% decrease in x^- Generally, one or two Gaussians proved 
sufficient to fit the broad Ha emission adequately. Figure [3] shows some examples of fits to 
Ha. Only a small number of galaxies with asymmetric profiles or very broad wings required 
a third Gaussian component. In most cases we were able to achieve an acceptable fit to the 
Ha+[N II] blend, but there are a few cases (~ 6) for which no acceptable fit to the narrow 
Ha line could be obtained with the [S II] model. For those cases, we relaxed the width 
constrains on the narrow lines. In each case, the width of the narrow Ha in the best-fit 
model remained quite close to the [S II] line width but the Ha fit was significantly improved 
by allowing its width to be a free parameter. 

We also modeled the H/3 + [O III] region in continuum-subtracted spectra, in the rest- 
frame range of 4800-5150 A, following procedures similar to those described above for the 
Ha-|-[N II] region. We used double-Gaussian models for the broad component of H/3 as well 
as for [O III], for which the double Gaussians represent a wing and a core component. The 
corresponding components in the two [O III] lines are constrained to have the same velocity 
width, and the wavelength separations and intensity ratios are fixed at their laboratory val- 
ues. For those galaxies with significant broad Fe II emission (mainly Fe II AA4924, 5018 AA), 



we ad opted an analytical model for the broad Lorentzian system "LI" in IVeron-Cetty et al. 



(l2004j . see their Appendix A). Each Fe II line was fitted with a Lorentzian profile, allow- 
ing the line center to vary within a small range around the expected value to account for 
shifting from the systemic redshift. The widths of both Fe II lines were constrained to be 
the same, while their fiux ratio was allowed to vary. In order to model the narrow B.f3, 



we follow iGreene &: Hoi (l2005bl ) to use the profile of [S II]. The H/3 centroid was fixed to 



the relative wavelength of narrow Ha if necessary, and i ts fiux was limited to be no larger 



than the value for Case B recombination (Ha = 3.1H/9, lOsterbrock &: Ferlandlll989l ). The 
procedure produced acceptable fit to H/3 in most cases, except for some objects (~ 13) in 
which the profile of [S II] seems not be a good model for narrow H/3, possibly resulting from 
small changes in spectral resolution across different echelle orders of the spectra. In these 
cases we used a single Gaussian to model narrow H/3, and we relaxed the width and fiux 
constraints. In Figure [3] we show some examples of best-fit models for the H/3 -|- [O III] 
region. Individual components for Fe II emission are also shown if they are present. It turns 
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out that the detected Fe II hnes in this sample are relatively narrow, typically 700 km s ^ 
FWHM, significantly narrower than the comm only used Fe II tenaplate I ZW 1, whose broad 



line system "LI" has FWHM = 1100 km s"^ flVeron-Cettv et all 120041 ) 



The FWHM of the overall profile of broad Ha is used for tracing the velocity dispersion 
of gas in the BLR, in the close environment of the black hole. We continue to use the 
FWHM(Ha) as a measurement of the line width for the whole sample in a consistent way, 
because the commonly used H/3 lines suffer from low S/N. Broad H/3 is often weak and 
sometimes not even detectable for objects in our sample. 

To estimate errors on the line width measurements, we should include statistical uncer- 
tainties from profile fits caused by noise and other random errors from spectral extraction 
and calibration, and uncertainties in continuum subtraction and deblending of broad com- 
ponents from narrow emission lines. In general, the formal fitting uncertainties from the 
profile fits are quite small, typically only 11 km s~^. Since the broad lines for this low- 
mass sample are narrow (mostly with FWHM < 2000 km s~^), the measurements of widths 
are much less aff ected by continuum subtraction than broad-lined AGN (see Section 2.5 of 



Dong et al.ll2008l ). Due to the small wavelength range over which the continuum-subtraction 
fits were performed and the good quality of the fits, the continuum subtraction does not 
add substantially to the error budget for the line widths. To explore the uncertainty from 
the emission- line profile debl ending in mor e detail, we create a set of artificial spectra from 



the best-fit model, following iGreene &: Hoi (l2005bf ). For each galaxy, we created a realiza- 



tion of the combined emission lines from the best-fit model parameters, and added Gaussian 
noise to match the S/N of the data. The artificial spectra created this way suffer from 
the same deblending difficulty as the original data. Then the artificial spectra are fitted 
with multiple-Gaussian models using the fitting procedure described above. The difference 
between the model and measured FWHM(Ha) is typically ~ 5% (or 0.02 dex). With these 
artificial spectra, we also investigated the uncertainty associated with the choice of model 
for the narrow-line profile, which often dominates the uncertainty in the decomposition of 
narrow/broad lines. We substituted the narrow-line model with the [S II] profile, the core 
of [O III], and narrow H/3 respectively, in the fitting procedure with the multiple-Gaussian 
model above. The typical standard deviation in the FWHM of broad Ha is ~ 8% (or 0.03 
dex). These two sources of error are combined to be typically ~ 9% (or 0.04 dex), and taken 
as the estimated uncertainty of FWHM(Ha). 
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3.3. Comparison of Broad-Line Widths With GH07 

It is instructive to compare our broad Ha measurements with those obtained from 
SDSS spectra by GH07. The smaller spectroscopic apertures for the Keck and Magellan 
data result in a smaller degree of starlight dilution of the AGN features, and the higher 
spectral resolution of the new data should permit more accurate deblending of the emission 
lines. Thus, we expect that the Keck and Magellan data should generally yield more accurate 
measurements of the broad Ha widths, particularly for very weak Ha emission lines. The 
SDSS spectra, on the other hand, have a more reliable flux calibration. In a small fraction 
of the objects from the GH07 sample, the broad Ha emission in the SDSS spectra is so weak 
that it is uncertain whether a distinct broad component is genuinely present or not, and our 
new spectra are particularly useful for testing the reality of these features tentatively seen 
in the SDSS data. 

We examined how our new measurements of the broad Ha FWHM compare with the 
results that GH07 obtained from fltting SDSS spectra, to test whether the objects with 
very weak broad Ha in GH07's low-mass AGN sample were genuine Seyfert 1 galaxies or 
not. We visually inspected the best-flt models for all the objects and divided them into two 
categories: objects with "deflnite" and "possible" broad Ha. In the following we refer to the 
two sub-categories as the d and p sub-samples, as indicated in Table [TJ For the d sample, 
the wings of broad Ha generally extend beyond the [N II] emission lines, and the broad Ha 
is signiflcantly wider than the narrow emission lines such as [N II]. We classifled objects as 
belonging to the p subsample if the peak amplitude of broad Ha was less than twice the RMS 
pixel-to-pixel deviation in the continuum-subtracted spectrum in the region surrounding Ha, 
or if the ratio of the flux of broad Ha to the RMS deviation of the continuum-subtracted 
spectrum was below 200. Some examples of "possible" broad Ha are shown in Figure HI 

Among the 93 galaxies in our sample, 14 (15%, 1 galaxy observed both by ESI and 
MagE) had ambiguous "possible" broad Ha emission. To investigate whether this was simply 
due to these 15 spectra having low S/N , we examined the average S/N in the Ha-|- [N II] 
region, and found that only three p objects had very low S/N of < 11. Overall, the objects 
classifled as being in the p category have a median S/N of 25, compared to a median value 
of 40 for our whole sample, so some of the p objects may simply be suffering from low S/N . 
However, some of the p objects have very high S/N spectra, so we can not attribute all the 
cases of ambiguity in broad Ha to low S/N . There are six p objects which only show very 
weak or possible broad Ha in the SDSS spectra (these are denoted as the c subsample in 
GH07). The other 8 out of the 14 p galaxies are classifled as broad- line AGNs by GH07. 
The difference between the SDSS results and our new fltting results could possibly be due 
to intrinsic AGN variability or variable AGN obscuration. Among the eight p galaxies. 
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SDSS J093147. 25+063503. 2 is confirmed to have outflow components in the narrow hnes, 
which might have mimicked part of the broad component in SDSS spectra. There are 
also five galaxies in the c ("candidate") sample of GH07 which we now classify as having 
definite broad Ha based on our new high-resolution spectra. A likely explanation is that 
the better quality of our new spectra makes the decompositions more accurate and better 
reveals the intrinsic properties of the broad components. Therefore, eight (9%) out of 82 
galaxies previously classified as broad-line AGNs by GH07 are found to have only possible 
broad Ha from the high-resolution spectra. This indicates the likely rate of false positive 
detections of broad Ha in the SDSS sample for objects near the threshold of detectability for 
broad Ha. In our sample, there is one object, SDSS J145045. 54— 014752.8, showing obvious 
double-peaked features in narrow lines, and also a definite broad component in Ha (see the 
last object of the right panel in Figure EJ. 

We compare our new measurements of broad Ha FHWM to those of GH07 by defining 
AFWHM = logFWHM(Ha)new - logFWHM(Ha)GH07- Figure E displays AFWHM vs. 
FWHM(Ha)ncw These two measurements are in reasonable agreement up to FWHM ~ 
1000 km s~^, but the difference increases as FWHM increases beyond 1000 km s~^. This is 
probably because the higher resolution of the new data enables us to fit the emission lines 
more accurately with more complex multi-Gaussian models, although some of the difference 
may be due to intrinsic source variability as well. 

4. Black Hole Masses and the Mbh — cr* Relation 

In this section we discuss our method of estimating black hole masses, and investigate 
the Mbh — cr^relation for our sample. In the following formulae for Mbh estimation, L5100 
denotes the AGN continuum luminosity XLx at A = 5100 A. 

4.1. Black Hole Mass 

The black hole mass is estimated via the virial relationship, Mbh = /-RblrAV^^/G, 
where -Rblr is the radius of the BLR, and the orbital velocity at that radius is estimated 
by the velocity width of the broad emission line, AV. Assuming that the gas in the BLR is 
virialized, the gas velocity traces the central mass in AGN, which is dominated by the black 
hole within the BLR radius. 



The broad line width can be measured with the FWHM or crune, whic h is the second 
moment of the profile, or the line dispersion. Both have merits and difficulties (jPeterson et al. 
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20041 ) ■ The line dispersion aunc has been suggested to be a more robust and precise estimator 
of viral velocity when measured from t he RMS spectra of reverberation mapping datasets 
(jPeterson et al.ll2004j : ICoUin et al.ll2006l ). However, it is highly sensitive to the contribution 



from the extended line wings, and therefore less rob ust in single-epoch s pectra when there is 



blending of other emission lines on the line wings (jPenney et al. 1120091). For measurements 



[1 



from single-epoch data, the FWHM is more commonly adopted (see a review in lMcGill et al. 



20081 ): it is sensitive to the line core and to the decomposition of the broad and narrow 
components, but it is relatively insensitive to the accuracy of measurement of faint extended 
wings on the broad-line profile. Denney et al. suggests that when the S/N is lower than 
10 — 20, both line width measurements would become unreliable, and line-profile fits would 
introduce systematic errors to single-epoch masses (e.g., ~ 0.17 dex offset in Mbh estimated 
with FWHM, their Table 5). In practice, we generally fit the broad Ha with two Gaussians. 
This is a purely empirical procedure and no specific physical meaning is assigned to the 
two components separately. Thus, we measure FWHM(Ha) from the overall profile of the 
Gaussian components used to fit the broad component of Ha. 



We follow GH07's approach (see iGreene fc Hoi l2007d . Appendix) to calculate Mbh 



which is estimat ed by using th e line-w idth of Ha and the BLR radius inferred from the broad 
Ha luminosity (IGreene &: Hdl2005bl ). In this work, the width of the broad Ha emission 
line is measured by decomposing the Ha+[N II] lines in the ESI and MagE data. The 
method relies on the br oad-line region radius-luminosity (-Rrlr — L) relation derived from 
reverberation mapping jKaspi et all boosi : ISentz et al.1 l20od . boool ) to determine the BLR 
radius from the estimat ed continuum l uminos ity, which in turn is estimated from the broad 
Ha luminosity following IGreene fc Hoi (l2007d ) since this is more accurately determined than 
the nonstellar continuum from our spectroscopic data. We update GH07's "recipe" fo r 
determining BH masses with the revised -Rrlr ~ L relation presented by lBentz et al.l (120091 ). 



who analyzed high-resolution HST images of 34 reverberation-mapped (RM) AGNs to obtain 
more accurate measurem ents of AGN continuum luminosity. The revised -Rblr ~ L relation 
from iBentz et al.1 koQ^ is 



log 



-Rblr \ 
It-days / 



i.5o;°:°^ + 0.5191°:°^^ log 



^5100 



10^4 ergs s" 



(2) 



We follow GH07 in assuming / = 0.75 ( lNetzerlll990l ) although we note that this / value 
is not derived for any specific physical model of the BLR. The choice of this part icular f factor 
aids in comparison of our results with prior work, as described below. From IGreene fc Ho 
( l2005bl ). the empirical relation between the line widths of broad Ha and H/3 is 
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3 /FWHM(Ha)^ i.03±0.03 



FWHM(H/3) = (1.07 ± 0.07) x 10^ ( - ----v--/ j 

\ lU Kill S / 

Combining these results with the virial relationship gives the BH mass as 



log 



M, 



BH 



Mo 



(4) 



If we substitute t he continuum l uminos ity with luminosity of Ha using the following 
empirical relation from I Greene fc Hoi (j2005bl ) , 



Lna = (5.25 ±0.02) x 10 



then we obtain the BH mass as 



42 



^5100 



1044 ergs s 



-1 



1.157±0.005 



ergs s 



(5) 



log 



M, 



BH 



Me 



6.40+°:°?+(0.45±0.05) log 



(6) 



Our Keck and Magellan spectra were not all taken under photometric conditions, and 
since the observations were obtained through narrow spectroscopic apertures, slit losses can 
be significant. The SDSS spectra have a more consistent flux calibration, so it is preferable 
to use L(Ha) measured from the SDSS data, even if this does introduce some additional 
uncertainty due to the fact that the Ha linewidths and luminosities are measured from non- 
simultaneous observations. Most objects in our sample are in the GH07 sample of active 
galaxies containing low-mass BHs, so we can obtain Lho for most of our sample from GH07. 
For the 5 obiects in BGH05 sam ple that were not included in GH07, we use the Ha luminosity 
Lhq fror nlGreene fc Hoi (120041) . One object in the BGH05 sample was not part of either the 
GH07 or lOreene fc h3 toO^ catalogs, and for this object we estimated Lho f rom the L.^mn 
measured by iBarth et al.l (120051 ) using the Lna-i^sioo relation (Equation [5]) of I Greene fc Ho 
(l2005bl ). 



The BH masses estimated by Equation El are systemat ically lower by 0.08 de x than those 



obtained using the method described in the Appendix of 
used the earlier version of the -Rblr — L relationship from 
to our mass estimator as Mghot; since it follows their basic method with only the -Rblr ~ L 



Greene fc Ho J2007ch . since they 
Bentz et al. We still refer 
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Greene fc Ho 



Greene &: Ho 



relatio nship updated. The RMS scatter in the Lua-L^ioo relation is ~ 0.2 dex ( 
2005bf) . If we take the uncertainties on the Ha luminosity as discussed in 
( l2005bl ). which are typically 0.13 dex, the typical formal uncertainties on the BH masses 
are 0.14 dex based on propagation of the measurement errors on the Ha luminosity and 
width. This random error does not include the important systematic uncertainty in the 
normalization factor /. This is comparab le to the observable e rror for BH masses in Seyfert 
galaxies (0.12 — 0.16 dex) estimated by iDenney et al.l ( l2009l ) based on single-epoch data, 
considering the effects of AGN variability and random measurement errors. 

In the discussion above, we followed GH07 and updated their recipe to calculate Mbh- 
The method is based on the Ha emission lines only (not requiring a continuum measurement) 
and it allows us to compare our results in a consistent way with previous work on SDSS AGNs 
presented by GH07. One particular point to consider is that our mass estimates assume a 
normalization factor / = 0.75 for the virial masses, for consistency with GH07 and other 
previous work, but this is not a unique choice for f. The virial factor has be en the subject 



of much discussion in the literature (e.g. lOnken et al.ll2004j : ICoUin et al. 



2006), and different 



/ values have been used in various recipes to obtain Mbh estimates. iMcGill et al.l (120081 ) 
compared Mbh estimators based on different emission-line and continuum measurements and 
showed that systematic errors could be as large as 0.38 ± 0.05 dex. 

We would like to examine how different recipes change the Mbh values of our sam- 
ple and the Mbh — o"* relation that will be discussed in Section 14.21 Most of the recipes 
make use of continuum luminosity, but we lack direct measurements of the optical AGN 
continuum luminosity for our sample. An alternative way is to substitute L.^mn with the 



broad Ha luminosity using the empirical relationship given by lGreene fc Hoi (l2005bl ). Some 
recipes rely on measurements of FWHM(H/3), and to test those relationships we substitute 
FWH M(H/3) with FWHM( Ha) using the empirical rel ationship between Ha an d H/3 w idths 
from I Greene fc Hoi ( l2005b[ ). If we use the recipe of IVestergaard fc Peterson! (12006) . the 
masses would be higher by about 0.25 dex. The formalism presented by IWang et al.l (|2009[ ) 
gives relatively larger BH mass than our estimator in the low mass end, typically about 0.6 
dex higher at BH mass of {log Mq^qy / Mq) = 6. They calibrated their BH mass estimator by 
fitting the reverberation-based masses for 35 AGNs using the continuum luminosity and H/3 
FWHM, as log(MBH/M0) = a-F7 log (Lgioo/lO^^ ergs s-i)+/3 log (FWHM(H/3)/103 km s-^) 
and fixed 7 to be 0.5. Their fits yield values of a = 1.39 and (3 = 1.09. This value of /3 is 
less than /3 = 2, which is commonly adopted for mass estimators. The 35 RM AGNs used 
in their fits mostly have BH masses in the range from 10^ Mq to 10^ Mq, and a tends to 
compensate for the mass difference in the fitting. This results in larger BH masses in the 
low mass range, i.e., for objects below about 10"^ Mq. 
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There is some lower limit to the BH masses that we are able to detect in SDSS data, 
due to a combination of factors. A possibly low BH occupation fraction in very low-mass 
galaxies would result in a low detection rate. The low luminosity of AGNs with small black 
holes makes the AGNs hard to detect. The host galaxy has to be bright enough to be 
spectroscopically targeted by SDSS, and the S/N of SDSS spectra is limited. Moreover, the 
large aperture of SDSS spectra can mix AGN emission with H II regions and dilute the AGN 
signal. Since we have identified the "definite" broad-lined AGNs that have genuine broad 
Ha in our sample, we obtain a cleaner sample that illustrates how low we can really go in 
selecting objects with low Mbh using SDSS. Our definite broad-lined sample with successful 
measurements of a* has a median BH mass (Mbh) = 9.5 x IO^Mq, and a minimum of 
2 X lO^Mp,. 



4.2. Mbh - (t* Relation 



Figure shows the Mbh ~ cr* relation for our sample of AGNs with low BH masses, 
and for active galaxies with higher black hole masses as well as nearby galaxies with direct 
dynamical measurements. The comparison sample of 56 active galaxies with higher BH 
masses based on sing le-epoch spectroscopy were selected from SDSS DR3 with z < 0.05 by 
Greene fc Hoi (l2006bl . hereafter GH06 sample), but the Mbh values have been re-calculated 
with our updated mass recipe (Equation E]). Literature data on 24 RM AGNs with stellar 
velocity dispersion measurements presented by IWoo et al.l (l2010l . and the references therein) 
were included. The Mb h for these RM AGNs were calculated from the virial products 
(VPs) listed in Table 2 of I Woo et al.l (120101 ). Since (7^ ^^ was used in VP, an isotropic velocity 
distribution gives f„ = 3, assuming aune = FWHM/2 (lOnken et al.ll2004l ): The ratio between 
FWHM and aune is different for different line profiles, i.e. FWH M/cr;j„s = 2.35 for a Gaussian 



profile, and it varies around an average of 2 (jCoUin et al.ll2006l ). Instead of the virial factor / 
obtained by Woo et al., we assume /a = 3 for consistency with the GH07 sample. This would 
decrease the masses by 0.24 dex compared to those listed in their work. We also included the 
well-kn own intermediate-mas s BHs in NGC 4395 with reverberation mass Mbh = (3.6±1.1) x 



lO^M^ (Peterson et al. 



from lFilippenko fc 



Ho rj2003h 



2005) and the velocity dispersion with an upper limit cr* < 30 km s ^ 
and P OX 52 with virial mass Mbh ^ (3.1-4.2) x lO^M© based 



on broad H/3 line width and L5100 (iBarth et al.l 120041: iThornton et al.l |2008| ) and velocity 



Barth et al. 



(I2004I). Masses for both objects wer e 



dispersion a* = 36 ± 5 km s^ from 
calculated utilizing the virial coefficient (/^r) = 5.5 ( lOnken et al.ll2004l : iPeterson et al.ll2004l ) 
we also scaled down both masses by 0.26 dex for consistency with other objects shown on 
the plot. Note that the virial normalization factor we assumed is arbitrary, and the derived 
Mbh was decreased by about 0.26 dex compared to masses derived using the / factor from 
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Onken et al.l (120041 ) . But the slope of Mbh — cr* relation is independent of /. Our low-mass 
sample appears to smoothly follow the extension of the Mbh— cr* relation for inactive galaxies, 
and there seems to be no strong change in slope across the whole mass range including active 
and inactive galaxies. We will quantify the slope for all the active galaxies described above. 

Assuming a log-linear form log Mbh = a + /31og(cr*/200 km s~^), we fit the slope and 
zero point of the Mbh — c"* relation for all the active galaxies with two regression meth- 
ods : the symmetric least- squares fitting method, fitexy (jPress et al.lll992l ) modified foUow- 
i nglTremaine et al.l (120021 ). and the maximum-likelihood estimate (MLE) method linmix_err 



(iKellyl 120071 ). both implemented in IDL. The former method accounts for uncertainties in 
both coordinates as well as the intrinsic scatter by adding a constant to the error in the 
dependent variable, and solves for the best linear fit by minimizing the reduced to unity. 
Linmix-err uses a Bayesian method to account for measurement errors and intrinsic scatter, 
and computes a posterior probability distribution function of parameters. Since there is no 
significant difference between the results from the two regression methods, we will only quote 
results from the MLE method linmix-err. We obtain a = 7.68 ± 0.08 and /3 = 3.32 ± 0.22, 
with eo = 0.46 ± 0.03 dex. The slope is a bit flatter than that for nearby ga l axies, w hich 
have = 4 . 24 ± 0.41 ( iGiiltekin et"aDl2009l ). but consistent with lOreene fc hJ (l2006bh and 
Woo et al.l ( 120101 ) who both showed evidence of a shallower slope for active galaxies than 
the inactive Mbh — o"* relation. If we fix the slope /3 to the best-fit value of 4.24 for nearby 
galaxies, we obtain, for the full sample of active galaxies, a zerop oint of a = 7.99 ± 0.0 4, 
which is -0.13 ± 0.09 offset from the value of a = 8.12 ± 0.08 from loiiltekin et al.l ( l2009l ). 



We now consider residuals in the Mbu — ct* relation, defined as AMbh = log(MBH/M0) — 
log ( Mbh /Mq) fit, where log(MBH/MQ)fit is calculated from a*, as a function of the bolo- 
metric luminosity and Eddington ratio (Figure [7]). We follow GH07's method for bolo- 
metric correction to estimate the bolometric luminosity from Lhq, according to Lboi = 
2.34 X 10^''(Lhq/10^^)°''^^ ergs s~^. Among the 93 objects in our sample, the median value of 
the Eddington ratio is (i^boi/-^Edd) = 0.3, where ^Edd = 1.26 x 1O^'^(Mbh/M0). The sample 
is dominated by objects radiating at substantial fractions of their Eddington limits, since 
the SDSS selection favors identification of the most luminous AGNs in any mass range. We 
find that the residual AMbh is significantly correlated with Lboi (Figure [7^). The Spearman 
rank correlation coefficient is = 0.33, with a probability P < 10~^ that no correlation 
is present. Moreover, AMbh shows a strong anti-correlation with Lboi/-^Edd (Figure [7)d, 
Vs = —0.53, P < 10^^). We are wary to overinterpret these correlations, since the Lboi/-^^Edd 
is formally anti-correlated with Mbh and Mbh is correlated with Lboi because both of them 
are deduced based on Lhq. 



Note that for our Mbh estimates we are using propagated measurement uncertainties 
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based on the errors in Ha width and luminosity, while the true uncertainties in Mbh are 
probably dominated by the uncertainty in the BLR geometry and the chosen value of /. 
If we adopt 3cr uncertainties in Mbh as measurement errors and do the regression again, 
there is little change in either the slope or zero point of the derived Mbh ~ cr* relation, 
but the intrinsic scatters decrease by about 35% to eo = 0.28 ± 0.05 dex. The regression 
parameters for di fferent samples and different uncertainties considered are listed in Table [2l 
Woo et al.l (120101 ) recently reported the Mbh ~ c"* relation for 24 RM active galaxies with 
BH mass 10^ < Mbh/M© < 10^. They obtained a slope /3 = 3.55 ±0.60 and intrinsic scatter 
eo = 0.43 ± 0.08, which are also consistent with our result. In the following, we will examine 
whether the intrinsic scatter varies as a function of mass. 

We divide the active galaxies mentioned above into two data sets, with lower and higher 
Mbh, and analyze the intrinsic scatter eo for the two data sets. The lower- Mbh data set 
includes our sarn ple and the Lick AGN Monitoring Project (LAMP) sample presented by 
Woo et al.l (I2OIOI ). The higher- Mbh data set in c ludes GH06 SDSS sample and 17 previous 
RM active galaxies (collected also in IWoo et al.l ( I2OIOI )). Assuming that the Mbh recipe we 
used is equally valid at all masses (which is not necessarily the case), we fit the two data 
sets and quantify the intrinsic scatter to be eo = 0.38 ± 0.04 dex for the Iow-Mbh sample 
presented here and eo = 0.40 ±0.04 dex for the higher- Mbh objects. The scatter for low-mass 
data set slightly changes when including NGC 4395 and POX 52, or excluding the LAMP 
sample. The scatter decreases a little to eo = 0.34 ± 0.04 dex if we exclude the "possible" 
broad Ha objects. The derived scatters for the two data sets are consistent, indicating that 
the intrinsic scatter in virial BH masses is not strongly mass- dependent. We also examine 
the standard deviation of residuals from the best-fit Mbh — o"* relation for the two data sets, 
and get 0.50 dex and 0.44 dex for the low and high mass range, respectively. 

The intrinsic scatter in Mbh ~ across the entire mass range is about 0.46 dex, cor- 
responding to a fa ctor of ~ 3. That is close to the intrinsic scatter for inactive galaxies, 
eo = 0.44 ± 0.06 (IGiiltekin et al.l boool ). We will show in gS] that the intrinsic scatter 
persists or even increases if we scale the active galaxies to follow the Mbh — c* relation of 
nearby galaxies with direct measurements of BH masses and then fit the relation for the 
combined sample of the active galaxies and nearby galaxies. Some of the dispersion in the 
Mbh ~ 0"* relation may result from different gal axy morphological types following different 
Ml 



cr^. loci, as recent ev i dence has suggested (|Hu 



2009; Giiltekin et al. 2009; Gadotti fc Kauffman 



3 



2oog 

2o"o9r 



I Greene et a. 



Greene et al 



2008; Graham fc Li 



2OIOI ). For example. 



the intrinsic scatter in Mbh ~ c"* relation of early-type g alaxies (referring to elliptical galaxie s 



and SO galaxies) is smaller than for late-type galaxi es (I Graham 



galaxy morphologies have been investigated by iGreene et al.l (120081 ) using HST images. The 



20081; 



Giiltekin et al. 



There is tentative evidence for a similar trend ii i the Greene fc Hoi (12004 ) sample, whose host 



20091). 
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difference could either be due to unaccounted systematic error in the Mbh measurements 
for spirals, or because the scatter in Mbh for spira ls is actually larger, as might happen if 
the behavior of BHs hosted by pseudobulges (see iKormendy fc KennicuttI 12004 for a re- 
view) and classical bulges are different. This has been suggested by IHuI ( 120081 ) . who found 
that pseudobulges tend to host lower Mbh than classical bulges at a given a*. Similarly, 
a recent study by iGreene et al.l ( l2010l ) presents stellar velocity dispersions for a sample of 
nine megamaser disk galaxies with accurately measured BH masses in the range of ~ 10^ 
Mq. They find that the na a ser ga laxies fall below the Mbh — cr* relation of elliptical galaxies 
defined by iGiiltekin et al.l ( 120091 ). Based on morphology and stellar population properties, 
they speculate that most of the nine megamaser galaxies are likely to contain pseudobulges, 
providing further support for the idea that the pseudobulge Mbh — tr* re lation is offset below 
the relation for classical bulges ( IHuI 120081 : iGadotti fc Kauffmannll2009l ). In §4.31 below, we 
investigate the possible dependence of the Mbh — c"* relation fit on host galaxy morphological 
type for our sample. 

Simulations of massive BH growth indicates that the scatter in the Mbh — cr* relation 
should increase toward the low-mass end. In models with high-mass seeds, galaxies evolve 
from well above the present-day Mbh — cr* relation toward the relation, while galaxies starting 
with low -mass seeds first lie far below the relation and their B Hs grow to move toward the 
relation (jVolonteril l2010l and references therein). Furthermore, IVolonteril ( 120101 ) suggested 
that for heavy seeds of ~ 10^ Mq, the low- mass end of the Mbh — cr* relation will flatten 
to an asymptotic value. In our sample, we do not find clear evidence for increasing scatter 
in the Mbh ~ o"* relation at low masses, and we can not discriminate among possible seed 
models based on our sample. Increasing scatter at the low-mass end should nevertheless be 
a generic result regardless o f the seed masses or the details of the BH growth mechanism. 
As discussed by iPengl (120071 ) . major mergers decrease the scatter of the Mnn-Mgai relation 
because of a central-limit tendency of the BH-to-galaxy mass relation to eventually converge 
to some mean value. 



Given the previous expectation that low-mass galaxies should show a larger scatter in 
the Mbh ~ c""* relation than high-mass galaxies, as well as the recent observational sugges- 
tions of an offset Mbh — o"* relation for pseudobulges relative to classical bulges, it remains 
somewhat puzzling that the Type 1 AGN samples do not show an obvious trend of either 
increasing scatter or changes in slope toward lower masses. Th is appears to be the case 
both for the small sample of AGNs with reverberation mapping ( IWoo et al.l 120101 ). and for 
the larger sample of objects with single-epoch masses presented here. It could be that the 
random errors on individual AGN virial BH mass estimates are so large that the scatter 
in the mass estimates (resulting from variations in BLR kinematics, inclination, or other 
properties) simply dominates over the intrinsic scatter in the Mbh — o"* relation for these ob- 
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jects, preventing us from detecting mass-dependent variations in the intrinsic scatter. Direct 
stellar-dynamical measurements of BH masses in nearby RM AGNs can be an important 
cross-check on the derived masses. However, while observations so far do not reveal evidence 
for any dramatic offsets between stellar-dynamical and reve rberation masses, only a few ob- 



jects are presently amenable to both types of measurements (jPavies et al.ll2006l : lOnken et al. 



2OO7I : iHicks fc Malk^boosh . For the vast majority of reverberation-mapped AGNs, the 



an- 



gular size of the BH's gravitational sphere of influence is too small to be resolved by either 
HST or by ground-based telescopes employing adaptive optics, and only a small number of 
broad-lined AGNs are currently suitable targets for stellar-dynamical measurements of Mbh- 



4.3. Morphology: Barred versus Unbarred Disk Galaxies 



Recent work by I Graham et al.l (1201 ll ) has suggested that the Mbh — cr* relation is differ- 
ent for AGNs with different host galaxy morphological types (i.e. barred disks or unbarred 
disks). For our sample, we can carry out a preliminary check for any offsets in the Mbh ~ cr* 
plane as a function of host morphology. The SDSS images are not sufficient to carry out 
a full census of host galaxy types for this sample, but there is morphological information 
for a subset of our sample from an HST W ide Field Planeta ry Camera2 (WFPC2) imag- 
ing snapshot survey with the F814W filter (j Jiang et al.ll201ll ). Among the galaxies in our 
sample having HST images, there are 38 objects which are morphologically classified as 
disk galaxies and which also have Mbh and velocity dispersion measurements, including 12 
barred and 26 unbarred galaxies. Only one galaxy from the GH06 higher-mass SDSS sample 
is included in the HST imaging survey; it is classifie d as an unbarred disk galaxy. Including 
the reverberation-mapped AGN sample compiled by lWoo et al.l ( l2010l ). as well as NGC 4395, 
and LAMP objects, we have morphological information for a total of 63 disk galaxy AGN 
hosts (25 barred and 38 unbarred). The morphological classifications for the reverberation- 
mapped ob j e cts ar e from NED where available, and otherwise adopted from clas sification by 
Bentz et al.l ( 120091 ). while the classification for LAMP objects are provided by iBentz et al. 
( 120111 . in preparation). If we remove the galaxies with a "possible" broad Ha component, 
there are 56 disk galaxies with morphological information (22 barred and 34 unbarred). One 
caveat in our division is that it is based on optical images, while it has been suggested that 
some galaxies appear to be unbarred in optical bands, but are very clearly ba rred when 



viewed in the infrared (e.g. iBlock fc WainscoatI Il991t iMulchaey fc ReganI 119971 ). Despite 



this caveat, we are limited by the data presently available and we carry out a preliminary 
examination based on the optical host-galaxy morphologies. 



We re-do the regression fits for the Mbh — o"* relation as in §4. 21 for different subsamples. 
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and the results are listed in Table |3l The best fit of the Mbh — cr* relation for the full 
subsample of disk galaxies (a = 7.50 ± 0.13, f3 = 3.04 ± 0.30) is consistent with that for 
the entire sample of active galaxies described previously (a = 7.68 ± 0.08, (3 = 3.32 ± 0.22). 
Within the uncertai nties, the slope i s cons istent with the slope of 4.05 ± 0.83 for 24 nonel- 
lipticals reported by iGiiltekin et al.l (120091 ). Fixing the slope to 4.05 yields an intercept of 
a = 7.89 ±0.06, which also agrees well with a = 8.01 ±0.16 for nonellipticals (IGiiltekin et al. 
2009[ l. The fit for barred disks only gives a = 7.81 ± 0.27 and /3 = 4.13 ± 0.72, similar to the 
Mny{ — relation (g = 7.80 ± 0.10, /3 = 4.34 ± 0.56) for barred inactive galaxies presented 
bv ICraham et al.l J201lh . 

Figure [8] shows the morphological type for the disk galaxies on the Mbh — cr* plot, and 
also the best fits of the Mbh — relations for barred and unbarred disks (see Table [3]). We 
consider our most reliable results to come from the subset of our sample that excludes the p 
objects with uncertain detections of broad Ha. For this subsample, we find a = 7.80 ± 0.29 
for barred disks versus 7.52±0.15 for unbarred disks, and the slope /3 = 4.03±0.80 for barred 
disks versus 3.01 ± 0.33 for the unbarred subset. Thus, the barred disks have a marginally 
larger zeropoint and steeper slope than unbarred disks, but the difference is not significant 
given the substantial uncertainties on th e fits. We also fit the M bh~o"^, relation with the slope 
fixed to that of inactive galaxies, 4.24 (IGiiltekin et al.ll2009l ). as an additional test for any 
offset. With fixed slope, the distinction between the barred and unbarred subsamples almost 
vanishes (see Table El) . Alternatively, if we were to adopt the high slope value proposed by 
Graham et al.l (120111 ) {(3 = 5.13), then the zeropoints would be all increased by about 0.35. 



4.4. Disk Inclination 

We have measured a* for our sample from the ESI and MagE spectra obtained with 
slit widths of 0'.'75 and 1", respectively. At the median redshift {{z) = 0.08, or luminosity 
distance Di=300 Mpc) for our sample, 1" corresponds to a scale of 1.5 kpc. If the bulge is 
not the dominant component of the host galaxy over that scale, then the slit spectra will 
be contaminated by the light from the disk, and the measured a* will include contributions 
from the rotational velocity of stars orbiting on the disk. Specifically, will be increased 
artificially in an edge-on disk, and might be slightly decreased in a face-on disk. In order to 
examine the effect of disk inclination on our sample, we use the axis ratio b/a to indicate 
inclination, where b/a close to 1 implies a face-on (low-inclination) system and b/a close to 
denotes an edge-on (high-inclination) system. The disk axis ratio is obtained from fits to 
HST images. For the subset of our sample having HST imaging, the axis r atio measurement s 



are based on the analysis of the WFPC2 snapshot imaging survey data (jjiang et al.ll201ll ). 
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For r everberation - mapp ed AGNs, axis ratios mea sured from HST im aginK data are taken 
from iBentz et al.l (120091 their Table 4) and from iBentz et al.l ( 120111 . in preparation). We 
divide the sample into three axis-ratio bins: b/a > 0.88, 0.72 < b/a < 0.88, and b/a < 0.72. 

Figure [9] shows the disk galaxies compiled in §4.31 in the Mbh — o"* plot, labeled by 
axis ratio. The low-inclination objects (large blue triangles) lie apparently to the left of the 
medium- (large orange squares) and high-inclination (large red circles) obj ects. We define 



AMrh as the Mbh deviation from the best-fit Mbh — o"* relation reported by iGiiltekin et al. 



( I2OO9I ). The distributions of AMbh for the three bins are shown in Figure [TUl A K-S test on 
the distributions of AMbh for the low- and high-inclination subsamples yields a probability of 
0.011, indicating that the distribution of the offset AMbh is likely to be significantly different 
for these two subsamples. Two extra PG objects (PG 1229+204 & PG2130+099) with Mbh 
and cr=K measurements are also shown on Figure M Their axis ratios of disk components are 
0.62 and 0.55, respectively, and they are classified as high-inclination objects. Including these 
two objects in the K-S test above does not alter the result. We also carry out fits to Mbh ~cr* 
relation for each subsample, and list the best-fit parameters in Table |H With the slope fixed 
to 4.24, the low-inclination subsample shows offsets from the other two subsamples by more 
than 0.4 dex. This gives additional evidence that the low-inclination and high-inclination 
systems are significantly offset in the Mbh — o"* plane, such that the more highly inclined host 
galaxies tend to have higher a^, at a given value of Mbh- The offset is consistent with the 
expectation that disk rotation contributes more significantly to the measured values of a* in 
the more highly inclined systems. This may be one source of error contributing to the scatter 
in the relation, and it may also introduce a bias to the zeropoint of the relation. Although 
it would be preferable to plot an Mbh — cr* relation that did not include this possible bias, 
there is no simple way to correct for the rotational contribution to the velocity dispersions 
from our data. 



4.5. Virial Normalization Factor 



In the discussions above, we have adopted an arbitrary virial factor / = 0.75 for con- 
sistency with most previous work on the SDSS-selected GH07 sample. Alternatively, we can 
also fit for the value of / that brings our sample into the best agreemen t wit h the local 
Mrh — cr* relation, following the same method used bv lOnken et al.l (|2004f) and lWoo et al. 



(120101 ) for reverberation-mapped AGNs. The values reported by IGiiltekin et al.l (120091 ) are 
a = 8.12 ± 0.08 and P = 4.24 ± 0.41. Fixing the slope to 4.24 for our full sample, we 
obtain a = 7.99 ± 0.04 for / = 0.75 as presented in §4.21 Therefore, the mean virial factor 
determined from this fit would be (/) = 1.01±0.21. Alternatively, fitting our sample to the 
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Mbh — c"* relation from iGraham et al.l ( 120111 ) would give a lower value of (/) = 0.52 ± 0.08. 
If we take the subsample of the 15 low-inclination galaxies above in §4.41 and scale to the 
Mbh — 0"* relation from Giiltekin et al., we get (/) = 0.56 ± 0.17. These results are reason- 
ably compatible with previous calibrations of (/) for reverberation-mapped AGNs, taking 
into account the fact that our work uses the FWHM as the measure of line width while 



Onken et al. 


2004; 


Woo et al. 


2010) 



relation is still limited by small-number statistics, and improving the determination of the 
virial normalization factor will require increasing the number of AGNs having high-quality 
reverberation mapping data as well as the number of spiral galaxies having direct measure- 
ments of Mbh from spatially resolved dynamics, and equally important, obtaining accurate 
a^: measurements for them. 

For completeness of our analysis, we fit the Mbh — o"* relation for the combined sample 
of active and nearby galaxies. First we use the virial factor derived above to recompute the 
BH masses for the active galaxies. Then we re-do the regression fit for the combined sample. 
Using the virial factor (/) = 1.01 derived above, which is based on normalizing the AGNs to 
the Mbh ~ relation of Giiltekin et al., we combine the AGNs with the 49 nearby galaxies 
with direct dynamical measurements of BH masses presented by the same authors. Then 
the regression for the combined sample yields a = 8.02 ± 0.05 and P = 3.81 ± 0.16 with an 
intrinsic scatter of eo = 0.46 ± 0.03. If the alternative (/) = 0.52 is used, then the fit yields 
a zero point of 7.89 ± 0.06 and a slope of 4.14 ± 0.16 with an intrinsic scatter of 0.49 ± 0.030 
The intrinsic scatter for the combined sample is the same as or even larger than that for 
active galaxies only. We note that the slope and the intrinsic scatter depend on how we 
select the virial factor, because the active galaxies are scaled to the local Mbh — o"* relation 
first and then included for the fit. 



^ Instead we cou ld choose (/} = 0.52 based on normalizing the AGNs to the Mbh — cr* relation of 
Graham et al.l (|201ll ). and then combine the AGNs with the 64 galaxies with direct BH mass measurements 
listed by Graham et al. Assigning a 10 percent error to the cr, for the 64 galaxies as Graham et al. suggested, 
the regression for the combined sample yields a = 7.93 ± 0.05 and (3 = 4.30 ± 0.16 with eo = 0.48 ± 0.03. 
Assigning a 5 percent error to the cr, instead will yield almost the same results. 
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5. Narrow Emission-line Properties 



5.1. Stellar vs. Gas Velocity Dispersion 

It is well-known that the velocity dispersion of the ionized gas in the narrow-line re- 
gion (NLR), most often measured from the [O III] A5007 A line, correlates with the stellar 
velocity dispersion of the host galaxy bulge, suggesting that gravity dominates the global 



kinematics of the NLR for active galaxies (e.g . IWhittle 



20031 : iGreene fc Hdl2005al : iBian et al.ll2006l : iKomossa fc Xu 



1985 



Nelson fc Whittle 
200?! : 



1996 



Barth et al 



Boroson 



20081 ). The 



high-resolution ESI and MagE data provides an opportunity to examine this relationship 
in galaxies of small a*. In addition to the [O III] line, we also consider lines from lower 
ionization species such as [N II] A6583 A and the [S II] AA6716, 6731 AA doublet, and check 
their relationship to cr* for our sample. 

For the [O III] AA4959, 5007 AA lines, each was fitted by a double- Gaussian model, 
representing a core and often (but not always) a blueshifted wing component. Generally 
the core component has a higher peak flux than the wing component. In some cases, the 
two-component model sometimes over-fits noise fluctuations near the line peak, resulting in 
individual Gaussian components with very narrow widths. We set an empirical criterion for 
the separation of the two components (> 80 km s~^) to define the blue- and red-shifts. The 
blue- or red-shifted wing component was removed from the FWHMcore measurements if the 
separation was > 80 km s~^. For the other cases, where the centroids of the two components 
were separated by less than 80 km s~^, we measured the FWHM from the [O III] profile 
and took the resulting value to be FWHMcore- Detailed studies of some galaxies where the 
NLR is spatially res olved reveal the possibility that the entire [O III] line may be outflowing 
(IFischer et al.l l201ll . and references therein). We use the relative velocity of the narrow 
component of H/3 to help determine the core component, since the Balmer recombination 
lines are less affected by the outflow when present. 

Figure [TT] illustrates the relationships between FWHM([0 III])/2.35 and ct=k for the core 
and wing components of the [O III] line as well as for the entire [O III] profile (Figure [Tib). 
The FWHM of the [O III] core correlates with cr*, with a Spearman rank correlation coeffi- 
cient Ts = 0.61, and a probability Pnuii less than 10~^ for the null hypothesis of no correlation 
(Figure [TTh). The relationship between the FWHM of the entire [O III] profile and cr» has 



a larger scatter, consiste nt with previous results (e.g. iNelson fc Whittle! Il996l : I Greene fc Ho 



2005al : iBarth et al.]l2008h . This illustrates that if we use the FWHM of [O III] as a proxy for 
cr* in active galaxies, the high-velocity wing component should be removed from the profile. 
The velocity shifts between the centroids of the wing to core components, defined as AVel, 
are mainly negative (Figure [T2l left panel), representing blue-shifts of the wing components. 
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while the centroid velocity of the core components are generally close to the systemic veloc- 
ity measured from the stellar ab sorption lines . If the blue-shifted win g component of [O III] 
is driven by winds or outflows (jWhittld Il985l : iKomossa et al.l 120081 ) . then the strength or 
velocity offset of the wing component could correlate with some measure of AGN activity, 
such as the Eddington luminosity ratio, Lboi/-Z^Edd- Thus we investigate how the equivalent 
width (EW) ratio of the [O III] wing to core components correlates with Eddington ratio, 
but we find no correlation between them. We find that AVel correlates with Lboi/-Z^Edd only 
weakly (r^ = —0.24, Figure [T21 middle panel). We also find that AVel correlates well with 
the FWHM of the wing component (r^ = 0.56, Pnuii < 0.001) and less strongly with the 
FWHM of the core component (r^ = OAl,Pnuu = 0.03). This confirms that the previously 
established correlation between [O III] blueshift and line width continues to hold in this 
low-mass regime. 

We also examine the correlation between gas and stellar velocity dispersions for the 
[S II] doublet and the [N II] A6583 A emission line, which has been decomposed from broad 
Ha. Since we use the [S II] profile as the model for the narrow lines for most cases (discussed 
in Sectioning), we can use either the [S II] or [N II] lines. The [N II] line is adopted here 
because there are a few cases for which [S II] is too weak to measure its width reliably. The 
FWHM is measured from the overall, combined profile of gaussians fitted to each narrow 
emission line. As shown in Figure [111 the FWHM of [N II] correlates well wit h cr » , with 
a Spearman rank correlation coefficient = 0.72 and Pnuii < 10~^. We follow iHol ( 120091 ) 



to define the gaseous velocity di spersion = FWHM/2.35 and the residual of ag — to 



be Act = logag — loga^,. While IHoI (|2009[ ) finds a correlation between Act and Eddington 
ratio for nearby active galaxies, we do not see such a correlatio n in our sa mple (Figure [121 
right panel) . Ou r sample most l y lies b elow the relation found by IHoI ( l2009l ) , but is generally 
consistent with [Greene fc Hoi (l2005al ). The objects scatter around Aa = (a^ = a^), 
regardless of increasing Eddington ratio. 



5.2. Emission-line Diagnostics 



We have measured emission lines of H/3, [O III] A5007 A, Ha, [N II] A6583 A, and [S II] 
AA6716, 6731 AA , which enable us to plot our sample on line-ratio diagnostic diagrams (e.g. 
BPT diagrams; Weilleux fc Osterbrockl [19871 ). as shown in Figure [T^ and to compare our 
measurements with the SDSS measurements obtained with larger spectroscopic apertures 
(fiber diameter of 3"). The flux of [O I] A6300 A was also measured from the starlight- 
subtracted spectra by fitting with a single Gaussian. From the new data (with aperture 
sizes of 0'.'75 x I'.'O and I'.'x (r.'5 ~ 3") for the ESI and MagE data, respectively), we see that 
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a few objects move out of the H II region of the BPT diagrams and into the Seyfert region 
above the maximum starburst hne ( iKewley et al.ll2006l ). With the small-aperture data, most 
objects lie in the Seyfert region of the diagram, with some objects extending to relatively 
low valu es of [N II] /Ha wh ich are indicative of metallicity lower than typical for classical 
Seyferts JCroves et aljl2006h . 



We investigate whether the line ratios measured from the small-aperture spectra sys- 
tematically change relative to measurements based on the SDSS spectra. The distributions 
of line ratios are shown on the top and right panels in Figure \572\ These distributions for the 
low-ionization ratios [N II]/Ha, [S Ilj/Ha, and [O I]/Ha only change slightly as a function of 
aperture size, with systematic offsets less than 0.08 dex between the median values for the 
small-aperture and SDSS measurements. This could be the effect of various parameters, due 
to the complex NLR conditions, including metallicity, ionization parameter, electron den- 
sity, and dust reddening. But the [O III]/H/3 ratio measured from the small-aperture data 
is systematically 0.3 dex higher than the values measured from the SDSS fiber aperture. 
The simplest interpretation is that in the larger SDSS aperture the emission lines included 
a larger contribution of emission from H II regions, which would have lower [O III]/H/3 than 
the AGNs. This would also imply that the host galaxies are actively star-forming, and we 
may be witnessing growth of the host galaxy bulge or pseudobulge coevally with the BH 
growth. Integral-field observations at high angular resolution would be particularly useful 
for examining the distribution and luminosity of H II regions in the host galaxies. 

These line-rat io diagrams are generally used to discriminate between star-forming galax - 
ies and AGNs (e.g. IHo et aLlEoOTl iKauffmann et al.lboosl : Iflao et al.l[2005[ iBarth et al.lboosh . 
However, as noted by GH07, the broad-line AGNs do not lie exclusively above the maximum 
starburst line of Kewley et al. on t he BP T diagrams, and some objects even fall below the 
empirical line of IKauffmann et al.l ( 120031 ) in the region of pure star-forming galaxies. That 
is, judging by the narrow-line ratios alone, these objects would not be classified as AGNs. 
For these objects, the detection of the broad Ha emission is the primary clue that an AGN 
is present, and if the broad emission lines were obscured by a dusty parsec-scale torus or by 
dust lanes in the host galaxy then these objects would not have been identified as AGNs at 
all. 

We explore whether the objects which fall below the maximum starburst line are unusual 
in terms of any of the measured parameters including a*, agas (indicated as FWHM([N Il])/2.35), 
the flux ratio between [O III] wing and core, the Fe II strength (relative to [O III]), 
FWHM(Ha), or the flux ratio between the narrow and broad components of Ha {Fn/b)- 
We find that this sub-population of objects does not exhibit any unusual properties in terms 
of these parameters, except in Fn/b-, which shows a higher relative strength of the narrow 
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components in this sub-population. Further investigation shows the median redshift of the 
sub-population {z ~ 0.1) is larger than the whole sample {z ~ 0.08). Thus the higher Fn/b 
ratio could be explained by the dilution of extended star formation contained in the aperture. 
It is unfortunately difficult to derive any detailed information on the stellar populations of 
the hosts from the spectra because of the strong AGN contribution to the spectra. We 
examined the ESI and MagE spectra, finding that higher-order Balmer absorption of signif- 
icant strength is visually apparent in 18 objects (about 20% of the sample), of which about 
10 objects are in the sub-population below the maximum starburst line. The absorption 
features imply contributions from A stars. Thus, compared with the sample as a whole, a 
larger fraction of objects in the H II sub-population appear to have significant contributions 
from intermediate-age stellar populations. 



6. Conclusions 

We have obtained new echelle spectroscopy of 76 Seyfert 1 galaxies selected to have low- 
mass black holes based on their broad-line widths and Lh^ luminosity, and we obtain reliable 
measuremen t s of cr ^, for 56 galaxies. Including the previous 17 Seyfert galaxies presented by 



Barth et al.l ( 120051 ). our low-mass AGN sample consists of 93 galaxies, of which 71 have 
measured velocity dispersions. The velocity dispersion ranges from 31 to 138 km s~^, and 
the data tend to lie on the extrapolation of the Mbh — cr* relation of inactive galaxies. We 
find that the intrinsic scatter in virial Mbh — o"* r elation is not strong ly mass-dependent. 



Combining our results with 56 SDSS active galaxies ( iGreene fc Hdl2006al ) and reverberation- 
mapped AGNs from the literature with measurements and Mbh estimation, as well as the 
two well-known intermediate-mass BHs in NGC 4395 and POX 52, we carry out new fits to 
the Mbh ~ c"* relation, finding a zero point of a = 7.68 ± 0.08 and slope (3 = 3.32 ± 0.22 for 
the Mbh — o"* relation in the form of log Mbh = a + /3 log(cr*/200km s~^), for an assumed 
virial normalization factor of / = 0.75. The AGN Mbh — cr* relation has an intrinsic scatter 
of 0.46 ± 0.03 dex, comparable to the intrinsic scatter of the Mbh ~ cr* relation observed for 
inactive galaxies. Among our low-mass BH sample with definite broad Ha components, we 
find that the lower limit of BH masses detectable in this SDSS sample is 2 x lO^M© (for 
/ = 0.75). We do not find a significant offset or slope difference in the Mbh — o"* relation 
between the subsamples of barred and unbarred disk galaxies. We do find that the disk 
galaxies with high inclination angles (edge-on systems) show a mild offset from the face-on 
systems in the Mbh ~ cr* relation. The rotation of the disk in the edge-on systems may 
artificially increase the measured a*, which may partly introduce a scatter in the relation. 
We also confirm that the narrow emission lines [N II] A6583 A, [S II] AA6716, 6731 AA, and 
the core of [O III] A5007 A (with the blue-shifted wing removed) have velocity dispersions 
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which trace the stellar velocity dispersion well, confirming that the forbidden emission- line 
widths can be used as a useful proxy for the stellar velocity dispersion in low-mass AGNs. 
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Table 1. Observations and Measurements 





z 


Exp. S/N 


fT(Mg lb) 


<T(Fe) 


fT(Call) 


IT* 




r VVxliVil nCK 


FWHM(Ha) log Lhoc 


log(MBH, 








(s) 




(km s ^) (km s 


^) (km s ^) (km s 


) 




(km s-i) 








SDSS JOOOlll. 15-100155.5 


0.0493 


5400 


38 


71 ± 3 






82+6 


76 + 


3 


1870 


1660 


40.11 


6.00 


m 


SDSS J002228. 36-005830. 6 


0.1059 


5400 


35 


54± 5 : 


60 + 


5 




57 + 


4 


807 


620 


41.37 


5.69 


m 


SDSS J004042. 10-110957. 6 


0.0277 


3600 


26 


57 ± 5 


55 + 


6 


55+5 


55 + 


3 


1530 


2380 


39.67 


6.13 


m 


SDSS J010712. 03+140844.9 ghOl 


0.0771 


2982 


16 


38 ± 4 


39 + 


5 


38+9 


38 + 


4 


914 


874 


41.44 


6.03 


b 


SDSS J011749. 81-100114. 5 


0.1411 


5400 


37 














756 


640 


41.53 


5.79 


m 


SDSS JO 12055. 92 -084945. 4 p 


0.1246 


5400 


14 


55 ± 6 


51 + 


6 




53 + 


4 


700 


1580 


40.82 


6.28 


m 


SDSS J014429.16-011047.3 


0.0609 


3600 


29 


68 ± 4 






73 + 12 


70 + 


6 


1760 


2520 


40.38 


6.50 


m 


SDSS J015804.75-005221.9 


0.0807 


3600 


34 


49+4 


40 + 


3 


47+8 


45 + 


3 


1680 


1180 


40.56 


5.90 


m 


SDSS .1022756.28+005733.0 


0.1280 


3600 


16 














929 


1040 


41.60 


6.25 


m 


SDSS J022849.51-090153.7 


0.0724 


5400 


43 


69 ± 3 






57 + 10 


63 + 


5 


697 


720 


40.39 


5.38 


m 


SDSS J023310.79-074813.3 p 


0.0312 


1800 


37 


101 ± 7 


114 + 10 


107+ 4 


107 + 


4 


1740 


1800 


39.92 


5.99 


m 


SDSS J024009.10+010334.5 


0.1956 


6300 


16 














737 


580 


41.64 


5.75 


m 


SDSS J024402.24-091540.9 


0.1220 


5400 


26 


73± 9 : 


79 + 


7 




76 + 


6 


970 


1000 


41.47 


6.16 


m 


SDSS J024912. 86-081525.6 gh02 


0.0297 


6000 


22 


49 ± 3 


52 + 


2 


60+7 


53 + 


3 


843 


702 


40.31 


5.32 


b 


SDSS .1030417.78+002827.3 


0.0450 


2400 


81 














1000 


900 


41.50 


6.08 


m 


SDSS .1032515.59+003408.4 gh03 


0.1023 


3000 


10 


54+6 


47 + 


7 




50 + 


5 


970 


886 


41.31 


5.98 


b 


SDSS J032707.32-075639.3 p 


0.1537 


3600 


22 


75+6 








75 + 


6 


697 


880 


40.80 


5.74 


m 


SDSS J034745.41+005737.2 


0.1792 


6600 


19 


107 ± 11 








107 + 


11 


865 


4717 


41.55 


7.58 


m 


SDSS J074836.80+182154.2 


0.0715 


5400 


28 


40+3 






45+8 


42 + 


4 


1660 


1760 


40.16 


6.08 


m 


SDSS J080629.80+241955.6 


0.0416 


1200 


28 


73+4 


71 + 


5 


69 + 12 


71 + 


5 


918 


1094 


40.83 


5.95 


e 


SDSS J080907.58+441641.4 


0.0541 


3600 


23 


66+4 


64 + 


3 




65 ± 


3 


950 


1104 


40.78 


5.94 


b 


SDSS J081550.23+250640.9 


0.0726 


3000 


17 


63+ 3 


67 + 


4 


65+5 


65 ± 


2 


903 


771 


40.63 


5.55 


e 


SDSS .1082325.91+065106.4 p 


0.0723 


3600 


9 


55+7 








55 ± 


7 


1320 


1300 


39.93 


5.70 


m 


SDSS J082347.95+060636.2 p 


0.1037 


3600 


13 


69+9 








69 + 


9 


1480 


1680 


40.21 


6.06 


m 


SDSS J082422.21+072550.4 p 


0.0815 


3600 


6 














1220 


960 


41.01 


5.92 


m 


SDSS J082443.28+295923.5 p 


0.0254 


1800 


47 


100+ 7 


104 + 


4 


118+ 6 


107 + 


3 


871 


691 


40.35 


5.33 


e 


SDSS J082912.67+500652.3 gh04 


0.0436 


2700 


35 


61+4 


62 + 


3 


58+4 


60 + 


2 


834 


759 


41.12 


5.76 


b 


SDSS .J083346. 04+062026.6 


0.1095 


5400 


26 


39+6 


51 + 


5 




45 + 


4 


1070 


640 


40.72 


5.42 


m 


SDSS ,1083928.45+082102.3 


0.1302 


3600 


38 














829 


600 


41.65 


5.78 


m 


SDSS .1084011.27+075915.7 p 


0.1324 


3600 


11 














1250 


1060 


41.09 


6.04 


m 


SDSS J090320.97+045738.0 


0.0567 


5400 


126 














784 


740 


41.50 


5.90 


m 


SDSS J090431. 21+075330.8 


0.0833 


5400 


89 














938 


860 


41.29 


5.94 


m 


SDSS .1091032.80+040832.4 p 


0.0732 


5400 


22 


72 + 12 








72 + 12 


864 


640 


40.09 


5.14 


m 


SDSS J091449.05+085321.1 


0.1398 


5400 


41 














849 


720 


41.68 


5.96 


m 



Table 1 — Continued 



idUoo iNiame 




z 


Exp. S/N 
(s) 


(T(Mg 

(km s~ 


lb) (T(Fe 
-1) (km s" 


tT(Ca II) (T, 
^) (km s~^) (km s~ 




FWHM(Ha) FWHM(Ha) log Lna 
GH07 (km s"!) 


log(MBH 


) KJOS 


SDSS J092547.32+050231.6 




0.1263 


5400 


22 


















760 


580 


41.55 


5.71 


m 


SDSS J092700. 53+084329.4 




0.1124 


5400 


26 


93 ± 12 : 


107 ± 15 






100 ± 10 


1150 


1220 


41.29 


6.26 


m 


SDSS J093147. 25+063503. 2 


p 


0.0853 


5400 


25 


52 ± 


9 










52 ± 


9 


755 


1460 


41.00 


6.29 


m 


SDSS J093147.25+063503.2 


p 


0.0857 


3000 


8 


40 ± 


8 


39 ± 


9 


28 ±12 


35 ± 


6 


755 


1360 


41.00 


6.22 


e 


SDSS J093829. 38+034826. 6 




0.1193 


5400 


28 


56 ± 


7 










56 ± 


7 


974 


800 


41.21 


5.84 


m 


SDSS J094057.19+032401.2 




0.0606 


5400 


210 


82 ± 


3 










82 ± 


3 


908 


800 


41.45 


5.95 


m 


SDSS J094310. 12+604559.1 


gh05 


0.0743 


3600 


16 


















807 


679 


41.34 


5.76 


b 


SDSS J094529. 36+093610. 4 




0.0131 


3600 


81 


75 ± 


2 


77 ± 


3 






76 ± 


2 


1930 


1720 


40.53 


6.22 


m 


SDSS J095151. 82+060143.7 




0.0932 


3600 


35 


70 ± 


6 


83 ±10 






76 ± 


6 


1260 


660 


41.01 


5.58 


m 


SDSS J100035.47+052428.5 




0.0785 


2100 


29 


















1000 


940 


41.62 


6.17 


m 


SDSS J101108.40+002908.7 


gh06 


0.1002 


1800 


10 


50 ± 


7 


61 ± 


7 






55 ± 


5 


1010 


1083 


41.36 


6.18 


b 


SDSS J101627.32-000714.5 


gh07 


0.0950 


3600 


15 






55 ± 


7 






55 ± 


7 




633 


41.06 


5.57 


b 


SDSS J102124.87+012720.3 




0.0668 


3600 


36 


78 ± 


3 










78 ± 


3 


1690 


1600 


40.39 


6.09 


m 


SDSS J102348. 44+040553.7 


p 


0.0988 


5400 


27 


91 ± 13 










91 ± 13 


696 


520 


40.95 


5.34 


m 


SDSS J103518. 74+073406. 2 




0.0674 


5400 


55 


109 ± 


4 










109 ± 


4 


867 


800 


41.27 


5.87 


m 


SDSS ,1104210.03-001814.7 




0.1144 


5400 


32 


















816 


800 


41.64 


6.04 


m 


SDSS J105755.66+482502.0 




0.0732 


1200 


22 


50 ± 


3 


49 ± 


4 


38 ± 


3 


45 ± 


2 


957 


898 


40.62 


5.68 


e 


SDSS J110501.97+594103.6 




0.0337 


900 


50 


124 ± 


5 


125 ± 


8 


120 ± 19 


123 ± 


7 


702 


3672 


40.51 


6.89 


e 


SDSS J110540.46+035309.0 




0.0993 


3600 


47 


















820 


800 


41.62 


6.03 


m 


SDSS J111031. 61+022043. 2 




0.0799 


5400 


92 


78 ± 


3 


76 ± 


4 






77 ± 


3 


1000 


920 


41.39 


6.05 


m 


SDSS J111749. 17+044315.5 




0.1082 


5400 


36 


62 ± 


9 : 


76 ± 


6 






69 ± 


5 


826 


680 


41.38 


5.77 


m 


SDSS J112526. 51+022039.0 




0.0490 


5400 


136 


91 ± 


2 






83 ± 10 


87 ± 


5 


1090 


940 


41.14 


5.96 


m 


SDSS J114339. 49-024316. 3 




0.0937 


5400 


61 


97 ± 


8 


98 ± 


5 






97 ± 


5 


919 


1460 


41.18 


6.37 


m 


SDSS J114343.76+550019.3 


P 


0.0272 


1800 


22 


31 ± 


3 


31 ± 


3 


32 ± 


3 


31 ± 


2 


1070 


1393 


40.39 


5.97 


e 


SDSS J114439.34+025506.5 


P 


0.1018 


5400 


21 


48 ± 


4 


47 ± 


6 






47 ± 


4 


942 


720 


40.51 


5.43 


m 


SDSS J114633.98+100244.9 




0.1245 


5400 


31 






62 ± 


8 






62 ± 


8 


790 


780 


41.54 


5.97 


m 


SDSS J115138.24+004946.4 


gh09 


0.1950 


1200 


13 


















810 


1304 


41.74 


6.52 


e 


SDSS J121518. 23+014751.1 




0.0713 


5400 


95 


74 ± 


4 






88 ± 


4 


81 ± 


3 


910 


1000 


41.08 


5.99 


rn 


SDSS ,1122.342.81+581446.1 




0.0146 


1800 


29 


44 ± 


3 


47 ± 


3 


46 ± 


3 


45 ± 


2 


979 


1577 


40.30 


6.04 


e 


SDSS ,1124035.81-002919.4 


ghlO 


0.0812 


1800 


26 


49 ± 


5 


63 ± 


6 


58 ± 


7 


56 ± 


3 


915 


713 


41.64 


5.93 


b 


SDSS J125055.28-015556.6 


ghll 


0.0815 


1800 


18 


68 ± 


5 


68 ± 


8 


64 ± 


8 


66 ± 


4 




2266 


41.27 


6.80 


b 


SDSS J131310.12+051942.1 




0.0492 


5400 


76 


68 ± 


3 






80 ± 


4 


74 ± 


3 


888 


580 


40.68 


5.32 


m 


SDSS J131310. 12+051942.1 




0.0489 


1200 


16 


64 ± 


5 


65 ± 


5 


68 ± 


6 


65 ± 


3 


888 


645 


40.68 


5.41 


c 


SDSS J131651.29+055646.9 




0.0554 


5400 


101 


82 ± 


3 










82 ± 


3 


1260 


980 


41.04 


5.95 


m 



Table 1 — Continued 





oUoo IN ame 


t lag 


z 


Exp. S/N 
(s) 


<j{Mg 
(km s 


lb) cr(Fe) (T(Call) cr, 

(km s"-*-) (km s~^) (km s~ 


) 


r VV rllVil^rl 


Cx.) r VV JrliVli^rlQ; j 


Og Lua 


10g(AJBH 


) KJDS. 


SDSS 


J131659. 37+035319. 8 


P 


0459 


5400 


36 


82 ± 


6 




80 ± 12 


81 + 


7 


887 


880 


40.61 


5.66 


m 


SDSS 


J131926. 52+105610.9 





0647 


5400 


174 


47 + 


3 






47 + 


3 


840 


860 


41.10 


5.86 


m 


SDSS 


J134144.51-005832.9 





1476 


5400 


31 














835 


660 


41.39 


5.75 


m 


SDSS 


J141234.67-003500.0 


ghl3 


1269 


1800 


16 














884 


945 


41.56 


6.15 


b 


SDSS 


J143450. 62+033842. 5 


ghl4 


0284 


1800 


17 


46 ± 


6 


63+3 


63+6 


57 + 


3 


1050 


1001 


40.34 


5.65 


b 


SDSS 


J144052. 60-023506.2 





0448 


5400 


143 








73+8 


73 ± 


8 


950 


840 


41.21 


5.89 


m 


SDSS 


J144705.46+003653.2 





0953 


5400 


44 


63 ± 


4 


65+7 




64 + 


4 


1160 


1500 


40.96 


6.29 


m 


SDSS 


J145045. 54-014752.8 





0996 


5400 


61 


131 ± 


6 


145 ± 10 




138 + 


6 


955 


2420 


41.49 


6.96 


m 


SDSS 


J150754.38+010816.7 





0613 


5400 


114 


131 ± 


5 




133+ 3 


132 + 


3 


699 


1680 


40.40 


6.14 


m 


SDSS 


J153425. 59+040806.7 





0395 


2200 


7 














927 


449 


40.02 


4.79 


e 


SDSS 


J154257.49+030653.2 





0655 


5400 


117 














1070 


980 


41.41 


6.12 


m 


SDSS 


J155005. 95+091035. 7 





0916 


5400 


51 


80 + 10 


77+8 




78 + 


6 


835 


1100 


41.10 


6.08 


m 


SDSS 


J161227.84+010159.7 





0973 


7200 


70 














944 


920 


41.40 


6.05 


m 


SDSS 


J161751. 98-001957.4 





0573 


5100 


14 


65 ± 


6 






65 + 


6 


1120 


1020 


40.37 


5.68 


m 


SDSS 


J162403. 63-005410.3 





0468 


5400 


35 


94 + 


4 






94 + 


4 


1150 


1340 


40.64 


6.05 


m 


SDSS 


J162636.40+350242.0 





0342 


1800 


34 


48 ± 


2 


53+2 


56+2 


52 + 


1 


802 


714 


40.76 


5.54 


e 


SDSS 


J163159. 59+243740.2 





0436 


2900 


38 


63 ± 


3 


61+3 


74+4 


66 + 


2 


839 


541 


40.84 


5.33 


e 


SDSS 


J170246. 09+602818.9 


ghl6 


0692 


1100 


18 






81 ± 11 


82+9 


81 + 


7 




1116 


41.25 


6.16 


b 


SDSS 


J172759. 15+542147.0 


ghl7 


0997 


3500 


17 






67+8 




67 + 


8 


806 


656 


41.27 


5.68 


b 


SDSS 


J205822. 14-065004.3 





0742 


5400 


62 


58 + 


3 : 






58 + 


3 


917 


860 


41.54 


6.06 


m 


SDSS 


J213728. 62-083823.3 





1609 


7200 


36 














904 


1000 


41.76 


6.29 


m 


SDSS 


J215658. 30+110343.1 





1081 


3600 


36 


88 + 


8 






88 + 


8 




990 


42.14 


6.45 


b 


SDSS 


J221139.16-010535.0 





0925 


7200 


20 


69 + 10 


68 + 11 




68 + 


7 


1070 


1540 


40.95 


6.31 


m 


SDSS 


J230649. 77+005023.4 





0610 


3600 


39 


58 ± 


3 




73+6 


65 + 


3 


1800 


1479 


40.55 


6.10 


m 


SDSS 


J232159. 06+000738.8 


ghl8 


1840 


3600 


12 


76 + 


9 






76 + 


9 




1531 


41.52 


6.56 


b 


SDSS 


J233837. 10-002810.3 


ghl9 


0357 


3600 


26 


58 ± 


3 


56+2 


55+3 


56 + 


2 




1553 


40.06 


5.92 


b 


SDSS 


J234807. 14-091202.6 





0779 


5400 


34 


80 + 


7 






80 + 


7 


1490 


1480 


40.76 


6.19 


m 



00 



Note. — Object names are given by their SDSS coordinate designations. Those names in iGreene fc Hoi b004) are also given. Objects with 
"possible" broad Ha are indicated as "p" . S/N is the mean signal-to-noise ratio per pixel in the spectrum in the spectral region around or 
redward of Mg Ift used to measure cr«. For MagE data, if the Mg lb fitting region and Fe fitting region are on the same echelle order, we 
measure a, from one region extending across these two, and list the result as o-(Mg lb). The ct* with ":" is dubious due to limi ted numbers 
of fitting templates. The FWHM(Ha)(3£fo7 ^Ha values are from GH07 if available and otherwise from I Greene fc Hoi ||2004| '1. For SDSS 



J215658. 30+110343.1, we use -L5100 from lBarth et all ||2005| ) and the Z/sioo-^Hq relation in lGreene fc Hoi ll2005bl) to e stimate Lng- Viria l mass 
estimates of Mbh are in units of Mq. The observational sub-samples are indicated as b, e and m for ESI sample in lBarth et al ] l|2005l 'l. new 
ESI observation and new MagE observation, respectively. 
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Fig. 1. — Comparison of MagE and ESI spectra (in black and red, respectively) for SDSS 
J093147.25+063503.2 and SDSS J131310. 12+051942.1. The ESI data have been rebinned 
to the spectral resolution of the MagE data. To facilitate comparison, all spectra have been 
rescaled to have flux density of unity at 5100 A. In each panel, the residual from subtraction 
of the MagE spectrum from the ESI spectrum (minus an arbitrary constant for clarity) is 
shown as a black line in the bottom. The blue horizontal line shows the arbitrary constant. 
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Table 2. Regression parameters 



Sample 


N 


Unccrt. 


a 


/3 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Full 


155 


la 


7.68 ± 0.08 


3.32 ± 0.22 


0.46 ± 0.03 


Pull 


155 


3a 


7.69 ± 0.08 


3.28 ± 0.22 


0.28 ±0.05 


Pull w/o "p" 


142 


la 


7.75 ± 0.08 


3.48 ± 0.21 


0.41 ±0.03 



Note. — Col. (1): Data set considered; the "Full" sample comprises 
objects with both Mbh and cr« available, including the SDSS sample and 
24 RM AGNs w ith stellar velocity dispersion measurements presented by 
IWoo et al.l \201(t }. NGC 4395 and POX 52. The "Full w/o p" subsample 
denotes the full sample minus the objects classified as having possible broad 
Ha. Col. (2): Number of objects in each sample. Col. (3): Adopted 
uncertainties for A/bh. Col. (4): Zero point assuming log Mbh = a + 
13 log((T,/200 km s"!), for / = 0.75. Col. (5): A/bh - cr, slope. Col. (6): 
intrinsic scatter. 



Table 3. Regression parameters for barred / unbarred disk galaxies 



Free fit Pit with fixed slope 



Sample 


N 


a 


/3 


eo 


a 


/3 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Disk 


63 


7.50 ±0.13 


3.04 ±0.30 


0.43 ± 0.05 


7.96 ±0.07 


4.24 


0.50 


Barred 


25 


7.81 ±0.27 


4.13 ± 0.72 


0.46 ± 0.09 


7.86 ± 0.09 


4.24 


0.40 


Unbarred 


38 


7.44 ±0.15 


2.79 ± 0.34 


0.42 ± 0.06 


8.04 ± 0.09 


4.24 


0.54 


Disk w/o "p" 


56 


7.56 ±0.13 


3.19 ±0.31 


0.42 ±0.05 


7.95 ±0.07 


4.24 


0.46 


Barred w/o "p" 


22 


7.80 ±0.29 


4.03 ± 0.80 


0.46 ± 0.09 


7.87 ±0.09 


4.24 


0.40 


Unbarred w/o "p" 


34 


7.52 ±0.15 


3.01 ± 0.33 


0.40 ± 0.06 


8.01 ±0.09 


4.24 


0.48 



Note. — Col. (1): Data set considered; the "Disk" ("Barred" / "Unbarred") sample comprises the 
"Full" sample with morphological type available and classified as (barred / unbarred) disk galaxies. The 
"Disk w/o p" subsample denotes the "Disk" sample minus the objects classified as having "possible" 
broad Hq. Col. (2): Number of objects in each subsample. Col. (3): Zeropoint assuming log A/bh = 
a +/3 log((7*/200 km s'^), for / = 0.75. Col. (4): A/bh - o"* slope. Col. (5): Intrinsic scatter. Col. 
(6): Zeropoint of A/bh ~ relation with the slope fixed to that of inactive galaxies listed in Col. (7), 
and the upper limits on intrinsic scatter listed in Col. (8). 
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Fig. 2. — Examples of ESI (left) and MagE (right) spectra of spectra in the Fe and CaT 
regions. The spectra are flux-calibrated in fx units and normalized to a flux level of unity. 
The observed spectrum is in black and the best-fitting broadened stellar template is in red. 



Table 4. Regression parameters for different disk inclination 



Free fit Fit witli fixed slope 



Sample axis ratio N a /3 eo a /3 eo 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



High b/a< 0.72 20 7.42 ± 0.25 3.24 ±0.70 0.50 ±0.10 7.74 ±0.11 4.24 0.48 

Medium 0.72 < 6/a < 0.88 20 7.70 ±0.20 3.86 ± 0.49 0.38 ±0.10 7.85 ± 0.09 4.24 0.36 

Low b/a> 0.88 15 7.59 ± 0.32 2.76 ± 0.67 0.44 ± 0.12 8.25 ±0.13 4.24 0.48 



Note. — Col. (1): Data set considered; the "High" ("Medium" / "Low") sample comprises the "Disk" sample 
with the axis ratio of the disk component available, divided into three bins of inclination angle. Col. (2): axis 
ratio range of each subsample. Col. (3): Number of objects in each subsample. Col. (4): Zeropoint assuming 
logMBH = a + 13 log{(Tt/200 km s^i), for / = 0.75. Col. (5): Mbh - O"* slope. Col. (6): Intrinsic scatter. 
Col. (7): Zeropoint of Mbh — relation with the slope fixed to that of inactive galaxies listed in Col. (8), and 
the upper limits on intrinsic scatter listed in Col. (9). 
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Fig. 3. — Examples of fitting to the H/3 + [O III] and Ha+[N II] regions for ESI (left) and 
MagE (right) data. The spectra are fiux-calibrated in fx units and arbitrarily normalized to 
a fiux level of unit around the peak of [O III] A5007 A. In each panel, continuum-subtracted 
spectra (black) and the best-fit model (blue) are shown, with their residuals shifted downward 
by an arbitrary constant for clarity. Individual components are overplotted: narrow lines 
(green), broad lines (magenta) and Fe II emission (cyan) if present. 



-42 - 



0.8 



0.4 



0.0 



1.0 



0.5 



0.0 



SDSS J023310.79 



I I Ill 



SDSS J1 14439.34 



SDSS J131659.37 



-4- 



0.4 



0.2 



0.0 



0.3 
0.2 
0.1 
0.0 



J- 












■ 1 1 








■ 1 

J 


' ' ' 1 1 1 1 1 1 


: A- -: 





4800 4900 5000 6500 6600 6700 



Rest Wovelength (A) 



Fig. 4. — Examples of galaxies showing possible broad Ha components. Fit components 
are illustrated as in Figure |3l 
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Fig. 5. — AFWHM(Ha), the logarithmic difference between the GH07 measurement of 
broad Ha FWHM and our new measurement, as a function of broad Ha hnewidth. Filled 
and open squares represent galaxies with definite and possible broad Ha, respectively. The 
inset box shows the distributions of broad Ha linewidths from GH07 and from our new 
measurements. 
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Fig. 6. — The Mbh — cr* relation for massive black holes. Filled dots represent new measure- 
ments from this work. Red dots are the MagE sample, purple dots are the newly observed 
ESI objects, and blue dots are the objects originally from the BGH05 sample. For objects 
observed with both ESI and MagE, both measurements are shown, connected by a short 
black line. The filled dots over-plotted with grey squares represent the galaxies with "pos- 
sible" broad Ha compo nents. Filled squares represent the 56 higher-mass active galaxies 
selected from SDSS by I Greene &: Hoi (l2006bl ). with BH masses updated with Equation [6l 



ments presented by 


Woo et al. 


( 


201 


o|). ^ 


Peterson et al. 


2005; 


Thornton et al. 


2008) 



NGC 4395 and POX 52 flFilippenko fc Hoi 12003 



to the virial factor of / = 0.75 or /o- = 3, in the virial mass Mbh = /-RblrAV^^/G (see 
text for deta ils). Open circ l es rep resent objects with dynamically determined BH masses 
compiled by iGiiltekin et al.l (120091 ) . and the red line is the Mbh — c* relat ion derived by 



Giiltekin et al 



J2009h . rled dashed-line is the Mbh — cr* relation reported by iGraham et al.l 



(120111 ). The blue solid line shows our best-fit Mbh ~ o"* relation for active galaxies with 
a = 7.68 ± 0.08 and /3 = 3.32 ± 0.22, and the blue dotted lines show the intrinsic scatter 
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Fig. 7.— The Mbh - a, residual, AMbh = log(MBH/M0) - log(MBH/M0)fit vs. (a) 
bolometric luminosity, (b) Lboi/-Z^Edd- Filled dots are the lo w-mass sample i n this w ork. Blue 
squares are those 56 active galaxies selected from SDSS by lGreene fc Hoi (l2006bl ). with BH 
masses updated with Equation [6l 
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Fig. 8. — The same as Figure [6l with the morphological classification of barred and unbarred 
disk galaxies shown. Small filled circles represent our sample with low-mass Mbh- Large 
magenta filled squares represent barred disk galaxies, while the cyan circles denote unbarred 
disks. Magenta and cyan lines show the best-fitting Mbh — cr* relation for the barred and 
unbarred subsamples that excludes the p objects with uncertain detections of broad Ha, 
respectively. 
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Fig. 9. — The same as Figure El enlarged to highlight the active galaxies classified as disk 
galaxies. Small filled circles, squares and crosses represent the same as Figure El Three 
subsamples of high-, medium- and low-inclination objects are overplotted with red open 
circles, orange open squares and blue triangles, respectively. The Mbh — cr* relation from 
Giiltekin et al. is shown for comparison. Only one galaxy (SDSS J0924+5607) from the GH06 
higher-mass SDSS sample is observed in the HST imaging survey. Two extra PG objects are 
also shown as stars. For PG 1229+204 and PG 2130+099, th e BH masses are 4 .0 x W Mq 
and 2.1 X IO^Mq (virial product from iPeterson et al.l l2004l : iGrier et al.l 120081 ) under the 
assumption of fa = 3 as discussed in §4.2[ and velocity dispersion are 162 ± 32 km s~^ and 



172 ± 46 km s ^ (jPasyra et al.l 120071 ). respectively. Their axis ratios of disk components are 
0.62 and 0.55, respectively, and they are classified as high-inclination objects. 
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Fig. 10. — Distribution of AMbh defined as the Mbh deviation from the best-fit Mbh — cr* 
relation from Giiltekin et al., in three bins of disk inchnation with the axis ratio range shown 
at the top right in each panel. From top to bottom is high- (close to edge-on), medium- 
and low-inclination (close to face-on). In each panel, the open and filled histograms show 
the distributions for all AGNs having host galaxy axis ratios (including the reverberation- 
mapped objects) and the subset from only our SDSS sample, respectively. 
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Fig. 11. — FWHM/2.35 vs. a^, for: (a) the core component of [O III] after removing the 
shifted wing component; (b) the blue-shifted or red-shifted wing component of [O III]; (c) 
the overall profile of [O III]; and (d) the overall profile of [N II]. The colors of plot symbols 
are as in Figure [61 The number at the top of each panel gives the Spearman rank correlation 
coefficient for the plotted data. The solid line in each panel represents FWHM/2.35 = cr,,. 
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Fig. 12. — Left panel: Distribution of velocity shifts relative to the stellar absorption system 
for the core component (black dashed) and the wing component (red dot-dashed) of [O III]; 
Distribution of AVel (km s~^), defined to be the relative velocity-shift between the wing 
to core components, is shown in blue solid histogram. Middle panel: AVel as a function 
of Eddington ratio. The object with |AVel| > 600 km s"^ is SDSS J145045.54-014752.8, 
showing obvious double-peaked features in narrow lines. Right panel: Distribution of the 
residuals Act = logcr^ — log ex* (cr^ = FWHM/2.35) versus Eddington ratio, Lf,oi/LEdd, for 
[N II] A6583 A. T he dashed line denotes cr^ = a^. The s o lid and dotted lines show the 
relations found by IHoI (l2009l . Equation 3) and lGreene fc Hoi (l2005al . Equation 4). 
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Fig. 13.— Diagnostic diagrams for [O III]/H/3 vs. [N II]/Ha, [S II]/Ha and [O I]/Ha. Filled 
circles represent objects with definite broad Ha, and open circles represent objects classified 
as having possible broad Ha. Grey triangles show the line ratios measured from the SDSS 
data by GH07 (open triangle s represent their ca ndidate sample). Dashed curves are the 
maximum starburst lines fromlKewley et al.l (120061 ): the solid curve is the pure star formation 
boundary ( Kauffmann et al. 20031 ): The horizontal dotted line marks [O III]/H/5= 3. 



